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Abstract 

A new class of phenomena discussed in this review is based on interaction between 
spatially separated, but closely located ferromagnets and superconductors. They are 
called Ferromagnet-Superconductor Hybrids (FSH). These systems include coupled 
smooth and textured Ferromagnetic and Superconducting films, magnetic dots, wires 
etc. The interaction may be provided by the magnetic fiux from magnetic textures and 
supercurrents. The magnetic fiux from magnetic textures or topological defects can pin 
vortices or create them, changing drastically the properties of the superconductor. On 
the other hand, the magnetic field from supercurrents (vortices) strongly interacts with 
the magnetic subsystem leading to formation of coupled magnetic-superconducting 
topological defects. We discuss possible experimental realization of the FSH. The pres- 
ence of ferromagnetic layer can change dramatically the properties of the supercon- 
ducting film due to proximity effect. We discuss experimental and theoretical studies 
of the proximity effect in the FSH including transition temperature, order parameter 
oscillations and triplet superconductivity. 
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1 Introduction 



In this review we discuss a new avenue in solid state physics: studies of physical phenomena 
which appear when two mutually exclusive states of matter, superconductivity and ferro- 
magnetism, are combined in an unified Ferromagnet- Superconductor Hybrid (FSH) system. 
In the hybrid systems fabricated from materials with different and even mutually exclusive 
properties, a strong mutual interaction between subsystems can dramatically change prop- 
erties of the constituent materials. This approach offers vast opportunities for science and 
technology. The interplay of superconductivity and ferromagnetism has been thoroughly 
studied experimentally and theoretically Pfl |2] for homogeneous systems. In such systems, 
both order parameters are homogeneous in space and suppress each other. As a result, one 
or both the orderings are weak. A natural way to avoid the mutual suppression of the order 
parameter of the superconducting (S) and ferromagnetic (F) subsystems is to separate them 
by a thin but impenetrable insulator film. In such systems the S and F subsystems interact 
via magnetic field induced by the nonuniform magnetization of the F textures penetrating 
into the superconductor. If this field is strong enough, it can generate vortices in the su- 
perconductor. The textures can be either artificial (dots, wires) or topological like Domain 
Walls (DW). The inverse effect is also important: the S currents generate magnetic field 
interacting with the magnetization in F subsystem. 

First experimental works on FSH were focused on pinning properties of magnetic dot 
arrays covered by a thin superconducting film IH 13 El EI ■ The effect of commensurability 
on the transport properties was reported in [SlIllElini- This effect is not specific for magnets 
interacting with superconductors and was first observed in textured superconducting films. 
First experiments with such films were performed in seventies. In these experiments the pe- 
riodicity of the vortex lattice fixed by external magnetic field competed with the periodicity 
of an artificial array created by experimenters. Martinoli et al. jHl El used grooves and 
Hebard et al. [HI [T2j used arrays of holes. This approach was further developed by exper- 
imentalists in nineties ^Bj-Pn]- Theoretical analysis was also performed in the last century 
|2U1 1^ 1^ . First observation of the dependence of vortex pinning by magnetic dots array 
on the magnetic field direction was presented by Morgan and Ketterson [7j. This was first 
direct indication of new physics in FSH. New insight into the FSH physics has been provided 
by Magnetic Force Microscope (MFM) and Scanning Hall Probe Microscope (SHPM). By 
using such imaging technique the group at the University of Leuven has elucidated several 
pinning mechanisms in FSH P^-|25j. 

Different mesoscopic magneto-superconducting systems were proposed and studied the- 
oretically: arrays of magnetic dots on the top of a SC film [23 HZl I2H1 1211; Ferromagnet- 
Super conductor Bilayer (FSB) |2H1 EOl ED E21 Ell , embedded magnetic nanowires com- 
bined with bulk superconductor 0^1 or superconductor film [271 [SHI , a layer of magnetic 
dipoles between two bulk superconductors [21], an array of magnetic dipoles mimicking the 
FM dots on SC film [10], "giant" magnetic dot which generates several vortices in bulk su- 
perconductor [H], single magnetic dots on a thin superconducting film [I^ li^ liH li^l iiO], 
thick magnetic film combined with thick [IT] UHl IHl ED] or thin superconducting film [^ I52j. 

The characteristic scale scale of the magnetic field and current variation in all mentioned 
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systems significantly exceeds the coherence length ^. It means that they can be considered 
in London approximation with good precision. In the next section we derive basic equations 
describing FSH. Starting from London-Maxwell equations, we derive a variational principle 
(energy) containing only the values inside either S or F components. These equations allowed 
us to study single magnetic dots coupled with superconducting film (Sec. 13. ip as well as 
arrays of such dots (Section 13. 2|) . The simplest possible FSH system - sandwich formed 
by Ferromagnetic and Superconducting layers, divided by ultrathin insulating film ( FSB),- 
can demonstrate unusual behavior: spontaneous formation of coupled system of vortices 
and magnetic domains. These phenomena are discussed in Section We also discuss the 
influence of the thick magnetic film on the bulk superconductor. 

An alternative approach to heterogeneous SC/FM systems is just to employ the proximity 
effects instead of avoiding them. The exchange field existing in the ferromagnet splits the 
Fermi spheres for up and down spins. Thus, the Cooper pair acquires a non-zero total 
momentum and its wave function oscillates in space. This effect first predicted by Larkin 
and Ovchinnikov jSHl and by Ferrel and Fulde [31] will be cited further as LOFF effect. 
One of its manifestation is the change of sign of the Cooper pair tunneling amplitude in 
space. At some conditions the Josephson current through a superconductor-ferromagnet- 
superconductor (S/F/S) junction has sign opposite to sirup, where (p is the phase difference 
between right and left superconducting layers. This type of junctions was first proposed 
theoretically long time ago by Bulaevsky et al. 110] and was called vr-junction in contrast 
to standard or 0-junction. It was first reliably realized in the experiment by Ryazanov and 
coworkers in 2001 |57| l58] and a little later by Kontos et al. jSHj- The experimental findings 
of these groups have generated an extended literature. A large exhausting review on this 
topic was published in the beginning of 2002 |60jj. A more special survey was published at 
the same time by GarifuUin ^1 . We are not going to repeat what was already done in this 
reviews and will focus presumably on works which appeared after its publication. Only basic 
notions and ideas necessary for understanding will be extracted from previous works. 

Most of the proximity phenomena predicted theoretically and found experimentally are 
based on the oscillatory behavior of the Cooper pair wave function. These are the oscil- 
lations of the transition temperature (first predicted in j62^. |63j), and the critical current 
vs. the thickness of ferromagnetic layer which are seen as oscillatory transitions from 0- to 
TT-junctions Other proximity effects besides the usual suppression of the order parame- 
ters include the preferential antiparallel orientation of the F-layers in a F/S/F trilayer, the 
so-called spin- valve effect jBH [HHl IHEI • 

More recently a new idea was proposed by Kadigrobov et al. [HZI and by Bergeret, Efetov 
and Volkov [68 : they have predicted that the magnetization varying its direction in space 
transforms singlet Cooper pairs into triplet ones. The triplet pairing is not suppressed by 
the exchange field and can propagate in the ferromagnet on large distances thus providing 
the long-range proximity between superconductors in S/F/S junctions. 

The proximity effects may have technological applications as elements of high-speed mag- 
netic electronics based on the spin valve action jSHI and also as elements of quantum com- 
puters inn]- Purely magnetic interaction between ferromagnetic and superconducting sub- 
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systems can also be used to design magnetic field controlled superconducting devices. A 
magnetic field controlled Josephson interferometer in a thin magnetic F/S bilayer has been 
demonstrated by Eom and Johnson [70] . 

In the next Section we derive basic equations. Third Section is focused on phenomena 
in FSH which are based on only magnetic interaction between ferromagnetic and supercon- 
ducting subsystem. Recent results on proximity based phenomena in bi- and tri-layer FSH 
are presented in the last Section. 

2 Basic Equations 

In the proposed and experimentally realized FSH a magnetic texture interacts with the super- 
current. First we assume that ferromagnetic and superconducting subsystems are separated 
by thin insulating layer which prevents proximity effect, focusing on magnetic interaction 
only. Inhomogeneous magnetization generates magnetic field outside the ferromagnets. This 
magnetic field generates screening currents in superconductors which, in turn, change the 
magnetic field. The problem must be solved self-consistently. The calculation of the vor- 
tex and magnetization arrangement for interacting, spatially separated superconductors and 
ferromagnets is based on the static London-Maxwell equations and corresponding energy. 
This description includes possible superconducting vortices. Londons approximation works 
satisfactory since the sizes of all structures in the problem exceed significantly the coherence 
length ^. We remind that in the Londons approximation the modulus of the order parameter 
is constant and the phase varies in space. Starting from the London-Maxwell equation in 
all the space, we eliminate the magnetic field outside their sources and obtain equations 
for the currents, magnetization and fields inside them. This is done in the subsection 2.1. 
In the subsection 2.2 we apply this method to the case of very thin coupled ferromagnetic 
and superconducting films. When proximity effects dominate, the Londons approximation 
is invalid. The basic equations for this case will be described in subsection 2.3. 

2.1 Three-Dimensional Systems. 

The total energy of a stationary F-S system reads: 



where B is the magnetic induction, M is the magnetization, Ug is the density of S-electrons, 
mg is their effective mass and v,, is their velocity. We assume that the SC density Ug and 
the magnetization M are separated in space. We assume also that the magnetic field B 
and its vector-potential A asymptotically turn to zero at infinity. Employing static Maxwell 
equation V x B = , and B = V x A, the magnetic field energy can be transformed as 
follows: 






(2) 
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Though the vector-potential enters exphcitly in the last equation, it is gauge invariant due 
to the current conservation divj = 0. When integrating by part, we neglected the surface 
term. This is correct if the field, vector-potential and the current decrease sufficiently fast at 
infinity. This condition is satisfied for simple examples considered in this article. The current 
j can be represented as a sum: j = js + jm of the SC and magnetic currents, respectively: 

nsfie 27r ^ , _ , 

js = i^V^-^A 3 
2ms $0 

j„ = cV X M. (4) 

We consider contributions from magnetic and S-currents into the integral (jSj) separately. We 
start with the integral: 

^Jj^AdV = iy'(V X M) ■ AdV (5) 
Integrating by part and neglecting the surface term again, we arrive at a following result: 

^j3^AdV = ^jM-BdV (6) 

We have omitted the integral over a remote surface /(n x M) ■ AdS. Such an omission is 
valid if the magnetization is confined to a limited volume. But for infinite magnetic systems 
it may be wrong even in simplest problems. We will discuss such a situation in the next 
section. 

Next we consider the contribution of the superconducting current to the integral (j21). In 
the gauge-invariant equation |21v2 is the phase of the S-carriers (Cooper pairs) wave-function 
and $0 = hc/2e is the flux quantum. Note that the phase gradient V(p can be included into 
A as a gauge transformation with exception of vortex lines, where (p is singular. We employ 
equation Q to express the vector-potential A in terms of the supercurrent and the phase 
gradient: 

A = — V(/7 ^j, 7 

271 HsC^ 

Plugging equation ((Tj) into equation we find: 

Since js = eTT-^v^, the last term in this equation is equal to the kinetic energy taken with the 
sign minus. It exactly compensates the kinetic energy in the initial expression for the energy 
Collecting all remaining terms, we obtain a following expression for the total energy: 

We remind again about a possible surface term for infinite magnetic systems. Note that 
integration in the expression for energy © proceeds over the volumes occupied either by 
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superconductors or by magnets. Equation (jH)) allows to separate the energy of vortices from 
the energy of magnetization induced currents and fields and their interaction energy. Indeed, 
as we noted earlier, the phase gradient can be ascribed to the contribution of vortex lines 
only. It is representable as a sum of independent integrals over different vortex lines. The 
vector-potential and the magnetic field can be represented as a sum of magnetization induced 
and vortex induced parts: A = Am, + A„, B = B^ + 'B^, where A^, (the index k is either 
m or v) are determined as solutions of the Londons-Maxwell equations: 

V X (V X Afe) = — j,, (10) 

c 

The effect of the screening of magnetic field generated by magnetization by superconduc- 
tor is included into the vector fields A^ and B^. Applying such a separation, we present 
the total energy Q as a sum of terms containing only vortex contributions, only magnetic 
contributions and the interaction terms. The purely magnetic part can be represented as 
a nonlocal quadratic form of the magnetization. The purely superconducting part is repre- 
sentable as a non-local double integral over the vortex lines. Finally, the interaction term 
is representable as a double integral which proceeds over the vortex lines and the volume 
occupied by the magnetization and is bi-linear in magnetization and vorticity. To avoid 
cumbersome formulas, we will not write these expressions explicitly. 



2.2 Two-Dimensional Systems. 

Below we perform a more explicit analysis for the case of two parallel films, one F, another 
S, both very thin and very close to each other. Neglecting their thickness, we assume that 
both films are located approximately at z = 0. In some cases we need a more accurate 
treatment. Then we introduce a small distance d between films which in the end will be put 
zero. Though the thickness of each film is assumed to be small, the 2-dimensional densities 
of S-carriers n^^^ = Ugds and magnetization m = M(i,„ remain finite. Here we introduced 
the thickness of the S film ds and the F film dm- The 3d super-carrier density ns(R) can 
be represented as ns(R) = 5{z)n'f'\v) and the 3d magnetization M(R) can be represented 
as M(R) = 5{z — (i)m(r), where r is the two-dimensional radius- vector and z-direction is 
perpendicular to the films. In what follows nf^ is assumed to be a constant and the index 
(2) is omitted. The energy © can be rewritten for this special case: 

where a = A(r, z = and b = B(r, z = 0). The vector-potential satisfies Maxwell-Londons 
equation: 

V X (V X A) = --A5{z) + -Vv5{z) (12) 

A mgC 

+ 47rV X {m5{z)) 
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Here A = \\/ds is the effective screening length for the S film, A^ is the London penetration 
depth and dg is the S-film thickness |7T|. 

According to our general arguments, the term proportional to Vv? in equation (|13|) de- 
scribes vortices. A plane vortex characterized by its vorticity q and by the position of its 
center on the plane Tq contributes a singular term to Vy?: 

V7 ( \ z^{v- rp) 

Vy^o r, Fq = q— ^ 13 

|r — ro| 

and generates a standard vortex vector-potential: 

g<l>o 5 X (r - ro) 



A„o(r - ro,z) 

ZTT I r — I'o I 
roo L (h\r - r„hp-^l^l 

dk (14) 



2n |r — rol 
°° Ji(A;|r - ro|)e-'=l^l 
1 + 2A;A 

Different vortices contribute independently into the vector-potential and magnetic field. A 
peculiarity of this problem is that the usually applied gauge divA = becomes singular 
in the limit dg, dm 0. Therefore, it is reasonable to apply another gauge Az = 0. The 
calculations are much simpler in Fourier-representation. Following the general procedure, 
we present the Fourier-transform of the vector-potential Ak as a sum Ak = A^k + A^,k. 
Equation for the magnetic part of the vector-potential reads: 

k(qA„k) - fc'A^k = ^ - 4mk x niqe*'^^^ (15) 

where q is projection of the wave vector k onto the plane of the films: k = kzZ + q. An 
arbitrary vector field Vk in the wave- vector space can be represented by its local coordinates: 

yk = V^z + VJlq + V^^izxq) (16) 
In terms of these coordinates the solution of equation ^T^i reads: 

= TT^e^'^' (17) 



1 ^ Am {kzm\ - qm^, 



ikzd 



Xk"^^ k^ ' ' 

Integration of the latter equation over kz allows to find the perpendicular component of 



_ AnXqim\ + tm^zj , 
"'-'I" l + 2Ag ' ^ ' 



whereas it follows from equation ()15|) that aj^^ = 0. Note that the parallel component of 
the vector-potential A^j^^^ does not know anything about the S film. It corresponds to the 
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magnetic field equal to zero outside the plane of F film. Therefore, it is inessential for our 
problem. 

The vortex part of the vector-potential also does not contain 2;-component since the 
supercurrents flow in the plane. The vortex solution in a thin film was first found by Pearl 
[72]. An explicit expression for the vortex-induced potential is: 

2i%{z X g)F(q) 

where F{q) = J2j ^^^''^ is the vortex form-factor; the index j labels the vortices and rj 
are coordinates of the vortex centers. The Fourier-transformation for the vortex-induced 
vector-potential at the surface of the SC film reads: 



i< !>o{z X g)F(q ) 
g(l + 2Ag) 



"^vq — ~ — , n \ V / 



The z-component of magnetic field induced by the Pearl vortex in real space is: 

271 Jo l + 2\q ^ ^ ^ ' 

Its asymptotic at z = and r A is B^z ~ $oA/ (vrr^); at r ^ A it is B^^ ~ *^'o/ (vrAr). Each 

Pearl vortex carries the flux quantum $o = vr^c/e. 

The energy (|TT| . can be expressed in terms of Fourier-transforms: 



H = + Hm + H^ra-, (23) 

where purely vortex energy is the same as it would be in the absence of the FM film: 
The purely magnetic energy H^a is: 

Hm = -lJ m^b^^ (25) 

It contains the screened magnetic field and therefore differs from its value in the absence of 
the SC film . Finally the interaction energy reads: 

Ughe f cPq 



-q(27r)2 
- ^/"^-^-(^ 

Note that the information on the vortex arrangement is contained in the form-factor F{q) 
only. 
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To illustrate how important can be the surface term, let consider a homogeneous perpen- 
dicularly magnetized magnetic film and one vortex in superconducting film. The authors 
PU] have shown that the energy of this system is = e^ — m^Q, where £° is the energy of the 
vortex in the absence of magnetic film, m is the magnetization per unit area and $0 = hc/2e 
is the magnetic flux quantum. Let consider how this result appears from the microscopic 
calculations. The vortex energy is just equal to e^. Purely magnetic term ()25|) does not 
change in the presence of vortex and is inessential. The first term in the interaction energy 
(I26p is equal to zero since the infinite magnetic film does not generate magnetic field outside. 
The second term of this energy is equal to — m$o/2. The second half of the interaction 
energy comes from the surface term. Indeed, it is equal to 

(l/2)hrny m{f x z) ■ Ard(p = —{l/2)j>A-dr 

= -m$o/2 

2.3 Eilenberger and Usadel Equations 

The essence of proximity phenomena is the change of the order parameter (Cooper pair wave 
function). Therefore, the London approximation is not valid in this case and equations for the 
order parameter must be solved. They are either Bogolyubov-DeGennes equations [ZSl 
for the coefficients u and v or more conveniently the Gor'kov equations [ZHI for Green func- 
tions. Unfortunately the solution of these equations is not an easy problem in the spatially 
inhomogeneous case combined with the scattering by impurities and/or irregular boundaries. 
This is a typical situation for the experiments with F/S proximity effects, since the layers are 
thin, the diffusion delivers atoms of one layer into another and the control of the structure 
and morphology is not so strict as for 3d single crystals. Sometimes experimenters delib- 
erately use amorphous alloys as magnetic layers [TH]. Fortunately, if the scale of variation 
for the order parameter is much larger than atomic, the semiclassical approximation can be 
applied. Equations for the superconducting order parameter in semiclassical approximation 
were derived long time ago by Eilenberger [77] and by Larkin and Ovchinnikov [TH]. They 
were further simplified in the case of strong elastic scattering (diffusion approximation) by 
Usadel [TS] . For the reader convenience and for the unification of notations we demonstrate 
them here referring the reader for derivation to original works or to the textbooks [SUl I81j. 

The Eilenberger equations are written for the electronic Green functions integrated in the 
momentum space over the momentum component perpendicular to the Fermi surface. Thus, 
they depend on a point of the Fermi-surface characterized by two momentum components, 
on the coordinates in real space and time. It is more convenient in thermodynamics to use 
their Fourier- components over the imaginary time, the so-called Matsubara representation 
[22] • The frequencies in this representation accept discrete real values Un = {2n + l)7rT, 
where T is the temperature. The case of singlet pairing is described by two Eilenberger 
anomalous Green functions F(co',k, r) and F^(u;,k, r) (integrated along the normal to the 
Fermi-surface Gor'kov anomalous functions), where uj stays for ujn-, k is the wave vector at 
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the Fermi sphere and r is the vector indicating at a point in real space (the coordinate of 
the Cooper pair center-of-mass) . The function F generally is complex in contrast to the 
integrated normal Green function (^(cjjk, r), which is real. Eilenberger has proved that the 
functions G and F are not independent: they obey the normalization condition: 

[G(^,k,r)p + |F(.;,k,r)p = l (27) 

Besides, the Eilenberger Green functions obey the following symmetry relations: 

F(cj, k, r) = F*{-u, k, r) = F\uj, -k, r) (28) 
G{-uj, k, r) = -G*(cj, k, r) = -G(cu, -k, r) (29) 



Eilenberger equations read: 

^ F(c<;, k, r) = 2AG(cj, k, r) 

j rfV(q)W^(k,q)[G(k)F(q) -F(k)G(q)] (30) 



( d 2e , , , 
2cj + V — -z— A r 
\ ar c 



where A(r) is the space (and time)- dependent order parameter (local energy gap); v is the 
velocity on the Fermi surface; Vr(k, q) is the probability of transition per unit time from 
the state with the momentum q to the state with the momentum k and p(q) is the angular 
dependence of the density of states normalized by / (fqp{(i) = N{0). Here iV(0) is the total 
density of states (DOS) in the normal state at the Fermi level. The Eilenberger equations 
have the structure of Boltzmann kinetic equation, but they also incorporate quantum coher- 
ence effects. They must be complemented by the self-consistency equation expressing local 
value of A(r) in terms of the anomalous Green function F: 



rji OO 

A(r)ln(-)+2vrT^ 



Air) 



UJr. 



- J d^kp{k)F{u;nXr) 



(31) 



In a frequently considered by theorists case of the isotropic scattering the collision integral 
in equation (|3U|) is remarkably simplified: 

J rfV(q)W^(k,q)[G(k)F(q) - F(k)G(q)] = ^ [G{k){F) - F{k){G)] , (32) 

where the relaxation time r is equal to inverse value of angular independent transition 
probability W and (...) means the angular average over the Fermi sphere. 
The Eilenberger equation is simpler than complete Gor'kov equations since it contains only 
one function depending on by one less number of arguments. It could be expected that 
in the limit of very short relaxation time Tcqt <^ 1 (Tco is transition temperature in the 
clean superconductor) the Eilenberger kinetic-like equation will become similar to diffusion 
equation. Such a diffusion-like equation was indeed derived by Usadel jT^j. In the case of 
strong elastic scattering and the isotropic Fermi surface (sphere) the Green function does 
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not depend on the direction on the Fermi sphere and depends only on frequency and the 
spatial coordinate r. The Usadel equation reads (we omit both arguments): 

2u}F - Dd (gSf - FdG) = 2AG (33) 

In this equation D = VpT/3 is the diffusion coefficient for electrons in the normal state and 
d stays for the gauge-invariant gradient: (9 = V — 2ieA/hc. The Usadel equations must be 
complemented by the same self-consistency equation (jHH). It is also useful to keep in mind 
expression for the current density in terms of the function F: 

j = ie2TTTN{0)D ^ {F*dF - FdF*). (34) 

a;„>0 

One can consider the set of Green functions G, F, F"^ as elements of the 2x2 matrix Green 
function g where the matrix indices can be identified with the particle and hole or Nambu 
channels. This formal trick becomes rather essential when the singlet and triplet pairing 
coexist and it is necessary to take in account the Nambu indices and spin indices simulta- 
neously. Eilenberger in his original article [77j has indicated a way to implement the spin 
degrees of freedom in his scheme. Below we demonstrate a convenient modification of this 
representation proposed by Bergeret et al. jH3I- Let us introduce a matrix ^(r,t; r',t') with 
matrix elements g^'^, , where n, n' are the Nambu indices and s, s' are the spin indices, defined 
as follows: 

g-Jiv^t-v'^t') = l^(f3)„,„. / rfe(V^n"s(r,t)V^L(r',t')) (35) 

n" 

The matrix fa in the definition is the Pauli matrix in the Nambu space. To clarify the 
Nambu indices we write explicitly what do they mean in terms of the electronic ■^-operators: 
ijjis = ijjg] ijj2s = "^1 ^iid s means —s. The most general matrix g can be expanded in the 
Nambu space into a linear combination of 4 independent matrices fk',k = 0,1,2,3, where 
tq is the unit matrix and three others are the standard Pauli matrices. Following [HS], we 
accept following notations for the components of this expansion, which are matrices in the 
spin space ^: 

9 = 90% + 93% + f-J = fiifi + /2if2 (36) 

The matrix / describes Cooper pairing since it contains only anti-diagonal matrices in the 
Nambu space. In turn the spin matrices /i and /2 can be expanded in the basis of spin Pauli 
matrices = 0, 1, 2, 3. Without loss of generality we can accept the following agreement 
about the scalar components of the spin-space expansion: 

/i = + /2<3-2; /2 = /oO"o + Ms- (37) 

It is easy to check that the amplitudes i = 0...3 are associated with the following combi- 
nations of the wave-function operators: 

^Each time when Nambu and spin matrices stay together we mean the direct product. 
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Thus, the amphtude /a corresponds to the singlet pairing, whereas three others are responsi- 
ble for the triplet pairing. In particular, in the absence of triplet pairing only the component 
/s survives and the matrix / is equal to 




Let us consider what modification must be introduced into the Eilenberger and Usadel equa- 
tions to take in account the exchange interaction of Cooper pairs with the magnetization 
in the ferromagnet. Neglecting the reciprocal effect of the Cooper pairs onto the electrons 
of d- or f-shell responsible for magnetization, we introduce the effective exchange field h{r) 
acting inside the ferromagnet. It produces pseudo-Zeeman splitting of the spin energy^. In 
the case of the singlet pairing the Matsubara frequency u must be substituted by u + ih{r). 
When the direction of magnetization changes in space generating triplet pairing, the Usadel 
equation is formulated in terms of the matrix g 

^d{gdg) - \uj\[f^a^,g] + sig-n.u[h,g] = -i[A,g], (38) 

where the operators of the magnetic field h and the energy gap A are defined as follows: 

/i = rga • h (39) 
A = Air2(T2 (40) 



To find a specific solution of the Eilenberger and Usadel equations proper boundary condi- 
tions should be formulated. For the Eilenberger equations the boundary conditions at an 
interface of two metals were derived by Zaitsev [HI] . They are most naturally formulated in 
terms of the antisymmetric [g"") and symmetric {g"^) parts of the matrix g with respect to 
reflection of momentum —pz assuming that z is normal to the interface. One of them 
states that the antisymmetric part is continuous at the interface {z = 0): 

= -0) = = +0) (41) 

The second equation connects the discontinuity of the symmetric part at the interface = 
g^{z = +0) — g^{z = —0) with the reflection coefficient R and transmission coefficient D of 
the interface and antisymmetric part g°- at the boundary: 

pmi - ) = - m 

^In reality the exchange energy has quite different origin than the Zeeman interaction, but at a fixed 
magnetization there is a formal similarity in the Hamiltonians. 
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where = g'^{z = +0) + g'^{z = —0). If the boundary is transparent (R=0, D=l), the 
symmetric part of the Green tensor g is also continuous. 

The boundary conditions for the Usadel equations, i.e. under the assumption that the 
mean free path of electron / is much shorter than the coherence length ^, were derived by 
Kupriyanov and Lukichev|HS]. The first of them ensures the continuity of the current flowing 
through the interface: 

. dg< . dg> 
"<^<^ = ">^>^' 

where the subscripts < and > relate to the left and right sides of the interface; (T<^> denote 
the conductivity of the proper metal. The second boundary condition connects the current 
with the discontinuity of the order parameter through the boundary and its transmission 
and reflection coefficients D{9) and R{9): 

, . dgy 3 ,cos6D(6).. , 

where 6 is the incidence angle of the electron at the interface and D{6), R{0) are correspond- 
ing transmission and reflection coefficients. This boundary condition can be rewritten in 
terms of measurable characteristics: 

where Rb is the resistance of the interface. In the case of high transparency {R. ^ 1) the 



boundary conditions ()43I44|) can be simplifled as follows 

df< dfy 



where 7 is the ratio of normal state resistivities. 



3 Hybrids Without Proximity Effect 
3.1 Magnetic Dots 

In this subsection we consider the ground state of a SC fllm with a circular very thin FM dot 
grown upon it. The magnetization will be considered to be flxed, homogeneous inside the 
dot and directed either perpendicular or parallel to the SC fllm (see flgureQ). This problem 
is basic one for a class of more complicated problems incorporating arrays of magnetic dots. 

We will analyze what are conditions for appearance of vortices in the ground state, where 
do they appear and what are magnetic flelds and currents in these states. The S fllm is 
assumed to be very thin, plane and inflnite in the lateral directions. Since the magnetization 
is conflned inside the flnite dot no difficulties with the surface integrals over infinitely remote 
surfaces or contours arise. 
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3.1.1 Magnetic Dot: Perpendicular magnetization 

For an infinitely thin circular magnetic dot of the radius R with 2d magnetization m(r) = 
mzQ{R — r)6{z — d) on the top of the SC film the magnetic field can be calculated using 
equations ()18ll9p . The Fourier-component of magnetization necessary for this calculation is: 

mk = z Ji{qR)e'^^'^, 47 

Q 

where Ji(x) is the Bessel function. The Fourier-transforms of the vector-potential reads: 

i^Ti'^mRJi^qR) 



A 



mk ^2 

X (e-i<^^^S^ + (e^fc.rf _ e-'?'^)) (48) 
^ l + 2gA ^ " ^ ' 

Though the difference in the round brackets in equation ()48p looks to be always small (we 
remind that d must be put zero in the final answer), we can not neglect it since it occurs 
to give a finite, not small contribution to the parallel component of the magnetic field 
between the two films. From equation pHjl we immediately find the Fourier-transforms of 
the magnetic field components: 

For the readers convenience we also present the Fourier-transform of the vector-potential at 
the superconductor surface: 

I iSn'^XmR ^ , , . 

In the last equation we have put e"^"^ equal to 1. 

Performing inverse Fourier-transformation, we find the magnetic field in real space: 



r. z) = 4.A™fl r Jj3mSI^,^i, (51) 

JO 1 + 2oA 



1 + 2gA 
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Figure 2: Magnetic field of dot witli and witliout vortex for R/X = 5 and ^o/STi'^'mR = 4 

B^^ (r, z) = 2nmR H MqR)Ji{qr)e-''\'\ 
Jo 

^ifh^^^'^ + " ~ (52) 

where Q(z) is the step function equal to +1 at positive z and —1 at negative z. Note that 
has discontinuities at 2; = and z = d due to surface currents in the S- and F-films, 
respectively, whereas the normal component B^ is continuous. 

A vortex, if appears, must be located at the center of the dot due to symmetry. If i? A, 
the direct calculation shows that the central position of the vortex provides minimal energy. 
For small radius of the dot the deviation of the vortex from the central position seems even 
less probable. We have checked numerically that central position is always energy favorable 
for one vortex. Note that this fact is not trivial since the magnetic field of the dot is 
stronger near its boundary. However, the gain of energy due to interaction of the magnetic 
field generated by the vortex with magnetization decreases when the vortex approaches the 
boundary. The normal magnetic field generated by the Pearl vortex is given by equation 
(jUJ. Numerical calculations based on equations (jSH 1221) ^^e case R> \ shows that B^ 
at the S-film {z = 0) changes sign at some r = Rq (see figure |2I) in the presence of the vortex 
centered at r = 0, but it is negative everywhere at r > i? in the absence of the vortex. 

The physical explanation of this fact is as follows. The dot itself is an ensemble of parallel 
magnetic dipoles. Each dipole generates magnetic field at the plane passing through the dot, 
which has the sign opposite to its dipolar moment. The fields from different dipoles compete 
at r < i?, but they have the same sign aX r > R. The SC current resists this tendency. The 
field generated by the vortex decays slower than the dipolar field vs. ). Thus, 

the sign of B^ is opposite to the magnetization at small values of r (but larger than R) 
and positive at large r. The measurement of magnetic field near the film may serve as a 
diagnostic tool to detect a S- vortex confined by the dot. To our knowledge, so far there were 
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Figure 3: Phase diagram of vortices induced by a magnetic dot. The hues correspond to 
the appearance of 1,2 and 3 vortices, respectively. 



no experimental measurements of this effect. 

In the presence of a vortex, energy of the system can be calculated using equations 
(I23ll26p . The appearance of the vortex changes energy by the amount: 

A = £^ + EjYiV (^3) 

where Sy = eoln^X/C,) is the energy of the vortex without magnetic dot, Eq = ^^/{IG'k'^X)] 
Emv is the energy of interaction between the vortex and the magnetic dot given by equation 
fr26|) . For this specific problem the direct substitution of the vector-potential, magnetic field 
and the phase gradient (see equations (j5(JI51|) )leads to a following result: 

^ f°° MqR)dq 
JO 1 + 2Xq 

The vortex appears when A turns into zero. This criterion determines a curve in the plane of 
two dimensional variables R/X and m^o/ey. This curve separating regimes with and without 
vortices is depicted in figure El The asymptotic of Emv for large and small values of R/X can 
be found analytically: 

R 

Emv ~ - m^o (y > 1) 
R R 

Emv ~ - "^*0^ (y < 1) 

Thus, asymptotically the curve A = turns into a horizontal straight line m$o/^t) = 1 at 
large R/X and logarithmically distorted hyperbola [m^^Q/ e.u){R/ X) = 2 at small ratio R/X. 

At further increasing of either itl^q/ev ot R/X the second vortex becomes energy fa- 
vorable. Due to symmetry the centers of the two vortices are located on the straight line 
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including the center of the dot at equal distances from it. The energy of the two- vortex 
configuration can be calculated by the same method. The curve 2 on figure (0) corresponds 
to this second phase transition. In principle there exists an infinite series of such transitions. 
However, here we limit ourselves with the first three since it is not quite clear what is the 
most energy favorable configuration for 4 vortices (for 3 it is the regular triangle). It is not 
yet studied what is the role of configurations with several vortices confined inside the dot 
region and antivortices outside. 

3.1.2 Magnetic Dot: Parallel Magnetization 

Next we consider an infinitely thin circular magnetic dot whose magnetization M is directed 
in the plane and is homogeneous inside the dot. An explicit analytical expression for M 
reads as follows: 

M = mo9{R- p)5{z)x (55) 

where R is the radius of the dot, mg is the magnetization per unit area and x is the unit 
vector along the x-axis. The Fourier transform of the magnetization is: 

Mk = 2'KmoR^^^^x (56) 

q 

The Fourier-representation for the vector-potential generated by the dot in the presence of 
magnetic film takes the form: 

= e^''&^UqR)cos{<P,) 

Let introduce a vortex-antivortex pair with the centers of the vortex and antivortex 
located at x = +po, x = —po, respectively. Employing equations (j23ll26|) to calculate the 
energy, we find: 

E = 2eoM-)-4eoAy^ TT^^' 

- 2rn.<,.Rr'-l^^^^^d, + E, (58) 
JO 1 + 2Xq 

where Eq is the dot self energy. Our numerical calculations indicate that the equilibrium 
value of Po is equal to R. The vortex-anti-vortex creation changes the energy of the system 
by: 

4 Jo 1 + ZXq 



- 2mo$oi? / — — TT dq (59) 

JO 1 + 2Ag 
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Figure 4: Phase diagram for vortices-anti-vortices induced by the magnetic dot with in-plane 
magnetization. 



The instabihty to the vortex-anti-vortex pair appearance develops when A changes sign. 
The curve that corresponds to A = is given by a following equation: 



The critical curve in the plane of two dimensional ratios and y is plotted numerically 
in figure (^. The creation of vortex-anti- vortex is energy unfavorable in the region below 
this curve and favorable above it. The phase diagram suggests that the smaller is the radius 
R of the dot the larger value is necessary to create the vortex-anti- vortex pair. At large 
values of R and mo$o > ^o, the vortex is separated by a large distance from the antivortex. 
Therefore, their energy is approximately equal to that of two free vortices. This positive 
energy is compensated by the attraction of the vortex and antivortex to the magnetic dot. 
The critical values of mo^o/^o seems to be numerically large even at R/X ~ 1. This is 
probably a consequence of comparably ineffective interaction of in-plane magnetization with 
the vortex. 

Magnetic dots with a finite thickness were considered by Milosevic et al. ^SIEHIE]- No 
qualitative changes of the phase diagram or magnetic fields were reported. 
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Figure 5: Left: Field dependence of the resistivity of a Nb thin film with a triangular array 
of Ni dots. (From Martin et al. P). 

Right: Critical current as a function of field for the high density triangular array at T = 
8.52 K Tc= 8.56 K (From Morgan and Ketterson 



3.2 Array of Magnetic Dots and Superconducting Film 
3.2.1 Vortex Pinning by Magnetic Dots 

Vortex pinning in superconductors is of the great practical importance. First time the 
artificial vortex pinning was created by S-film thickness modulation in seventies. Martinoli 
et al. [inj have used grooves on the film surface to pin vortices and Hebard et al. [TH [T2] 
have used triangular arrays of holes. Magnetic structures provide additional possibilities to 
pin vortices. First experiments were performed in the Louis Neel lab in Grenoble [3, 4J. 
These experiments were performed with dots several microns wide with the magnetization 
parallel to the superconducting film. They observed oscillations of the magnetization vs 
magnetic field. These oscillation was attributed to a simple matching effect: pinning becomes 
stronger when vortex lattice is commensurate with the lattice of pinning centers. This can be 
measured in terms of external, normal to the film magnetic field needed to generate integer 
number of vortices per unit cell of the pinning array. 

Flux pinning by a triangular array of submicron size dots with typical spacing 400-600nm 
and diameters close to 200nm magnetized in-plane was first reported by Martin et al. [6]. 
Oscillations of the resistivity with increasing flux were observed with period corresponding to 
one flux quanta per unit cell of magnetic dot lattice (see figure EJjeft. This can be explained 
by the matching effect. Though matching effect is not specific to magnetic pinning arrays, 
enhanced pinning with magnetic dots with magnetization parallel to the film was observed 
by Martin et al. jUj. 

Dots array with out-of-plane magnetization component was first created and studied 
by Morgan and Ketterson f7j. They have measured critical current as a function of the 
external magnetic field and found strong asymmetry of the pinning properties vs magnetic 
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Figure 6: Schematic presentation of the polarity dependent flux pinning, presenting the 
cross section of a Pb film deposited over a magnetic dipole with in-plane magnetization: (a) 
A positive FL (wide gray arrow) is attached to the dot at the pole where a negative flux 
quantum is induced by the stray field (black arrows), and (b) a negative FL is pinned at the 
pole where a positive flux quantum is induced by the stray field. (From Van Bael et al. J87|) 



field direction (see figure ERight). This experiment has given the first direct experimental 
evidence that the physics of vortex pinning by magnetic dots is different from that of common 
pinning centers. 

Pinning properties of the magnetic dots array depends on several factors: magnetic 
moment orientation, the strength of the stray field, the ratio of the dot size and the dot lattice 
constant to the effective penetration depth, array magnetization, the strength and direction 
of the external field, etc. The use of magnetic imaging technique, namely Scanning Hall 
Probe Microscope (SHPM) and Magnetic Force Microscope (MFM) has revealed exciting 
pictures of vortex "world" . Such studies in combination with traditional measurements gives 
new insight in vortex physics. This work was done mainly by the group at the University of 
Leuven. Below we briefly discuss only a few cases studied in great details by this group. 

Dots with Parallel Magnetization. Van Bael et a/.j^ studied with Scanning Hall Probe 
Microscope (SHPM) the magnetization and vortex distribution in a square array (1.5/im pe- 
riod) of rectangular (540nmX360nm) cobalt trilayer Au(7.5nm)/Co(20nm)/Au(7.5nm) dots 
with magnetization along the edges of the dots lattice. SHPM images have revealed mag- 
netic field redistribution due to superconducting transition in the covered 50nm thin lead 
superconducting film. These data were interpreted by Van Bael et a/.jHZI as formation of 
vortices of opposite sign on both sides of the dot. By applying external magnetic field Van 
Bael et a/.jSTj have demonstrated the commensurate lattice of vortices residing on the "end" 
of magnetized dots. This location is in agreement with theoretical prediction [351. Remark- 
ably, they were able to observe "compensation" of the vortices created by the dots stray field 
with vortices of the opposite sign due to applied normal field (see figure IH)). 

Dots with Normal Magnetization. Van Bael et a/.jSH] have elucidated with the SHPM 
images the nature of previously reported (see e.g. work [7j) anisotropy in the vortices pinning 
by the array of dots with normal magnetization. They have used 1/xm period lattice of square, 
400nm side length and 14nm thin, Co/Pt multilayer dots covered with 50nm thin lead film. 
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Figure 7: SHPM images of a (10.5/im)^ area of the sample in H=-1.6 Oe (left panel) and 
H=1.6 Oe (right panel), at T=6.8 K (field-cooled). The tiny black/white dots indicate the 
posi- tions of the Co/Pt dots, which are all aligned in the negative sense (m < 0). The flux 
lines emerge as diffuse dark (H < 0) or bright (H > 0) spots in the SHPM images. (From 
Van Bael et a/.|88i) 



Zero field SHPM images show the checkerboard-like distribution of magnetic field (see Sec. 
3.3.4) The stray field from the dots were not sufficient to create vortices. In a very weak (1.6 
Oe) external field the average distance between vortices was about 4 lattice spacings. In the 
case of the field parallel to the dots magnetization vortices reside on the dots, as the SHPM 
image shows (see figure [7^). In the case of the same field with opposite direction, the SHPM 
shows vortices located at interstitial positions in the magnetic dots lattice (see figure [7|d). It 
is plausible that the pinning barriers are lower in the second case. 

Figure |H1 shows dependance of superconduction film magnetization versus applied mag- 
netic field normal to the film. Moshchalkov et al. jSl] have shown that magnetic field depen- 
dence of film magnetization of the superconducting film is very similar to the critical current 
dependence on magnetic field. Figure |H1 shows strong anisotropy of the pinning properties on 
the external magnetic field direction. Magnetic field parallel to the dots magnetic moment 
shows much stronger vortex pinning than antiparallel. 

3.2.2 Magnetic Field Induced Superconductivity 

Consider a regular array of magnetic dots placed upon a superconducting film with mag- 
netization normal to the film. For simplicity we consider very thin magnetic dots. Namely 
this situation is realized magnetic films with normal magnetization used in experiment |90j ) . 
The net flux from the magnetic dot through any plane including the surface of the supercon- 
ducting film (see figure Ej) is exactly zero. Suppose that on the top of the magnetic dot the 
z-component of the magnetic field is positive as shown in the mentioned figure. Due to the 
requirement of zero net flux the z-component of the magnetic field between the dots must 
be negative. Thus, S-film occurs in a negative magnetic field normal to the film. It can be 
partly or fully compensated by an external magnetic field parallel to the dot magnetization 
(see figure ini). Such a compensation can be even more effective in for a regular array of 
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Figure 8: M (H/Hi) magnetization curves at different temperatures near Tc (7.00 K open 
symbols, 7.10 K filled symbols) showing the superconducting response of the Pb layer on top 
of the Co/Pt dot array with all dots aligned in a positive upper panel and negative lower 
panel sense. Hi=20.68 Oe is the first matching field. (From Van Bael et a/.|88^) 
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Figure 9: Schematical magnetic field distribution in the the array of dots with normal to 
the superconducting film magnetization. (From Lange et al. cond-mat/0209101j) 

magnetic wires embedded in alumina template [SHI IHZl EI] • Lange et al. 100] have proposed 
this trick and reached a positive shift of the S-transition temperature in an external mag- 
netic field, the result looking counterintuitive if one forgets about the field generated by the 
dots. In this experiment a thin superconducting film was covered with a square array of 
the CoPd magnetic dots with normal to the film magnetization. The dots had square shape 
with the side 0.8/im, the thickness 22nm and the dot array period 1.5/im. The H-T phase 
diagrams presented in [90J for zero and finite dots magnetization demonstrate appearance 
of the superconductivity by applying magnetic field parallel to the dot magnetization. At 
T=7.20K the system with magnetized dots is in normal state. It undergoes transition to the 
superconducting state in the field 0.6mT and back to the normal state at 3.3mT. From the 
data in figure 3 in work by Lange et a/.|H3 one can conclude that the compensating field is 
about 2mT. 

3.2.3 Magnetization Controlled Superconductivity 

Above (Sec. I3.2.2j) we have discussed example when application of magnetic field can trans- 
form FSH system from normal to superconducting state. This was due compensation of the 
dots stray magnetic field with external magnetic field. 

Earlier Lyuksyutov and Pokrovsky [2E] have discussed theoretically situation when de- 
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/MAGNETIC DOT SUPERCURRENT \ /MAGNETIC DOT SUPERCURREN' 

Figure 10: Left: Magnetized magnetic dots array. Vortices of different signs are shown 
schematically by supercurrent direction (dashed lines). The magnetic moment direction 
is indicated by ±. Both vortices bound by dots and created spontaneously are shown. 
Magnetized array of dots create regular lattice of vortices and antivortices and provide strong 
pinning. Right: Demagnetized magnetic dots array results in strongly fluctuating random 
potential which creates unbound antivortices/ vortices, thus transforming superconducting 
film into resistive state. 
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Figure 11: Left:The checker-board average structure of the vortex plasma. Right: The 
average number of the unbound vortices in the cell of size a via the parameter k propor- 
tional to the dot magnetic moment. Dot-dashed line corresponds to T/eo = 0.15, solid line 
corresponds to T/eo = 0.4, dashed line corresponds to T/eo = 2. 

magnetized array of magnetic dots with normal magnetization create resistive state in the 
coupled superconducting film. However, superconducting state can be restored by magneti- 
zation of the dots array. This counter intuitive phenomena can be explained on qualitative 
level. In the case when single dot creates one vortex, magnetized array of dots results in 
periodical vortex/antivortex structure with anti-vortices localized at the centers of the unit 
cells of the square lattice of dots as shown in figure ITULeft. Such order provides strong pin- 
ning. More interesting is demagnetized state in which the induced vortices and antivortices 
create a random field for a probe vortex. If the lattice constant of the array a is less than the 
effective penetration depth A, the random fields from vortices are logarithmic. The effective 
number of random logarithmic potentials acting on a probe vortex is = (A/a)^ and the 
effective depth of potential well for a vortex (antivortex) is y/Ney. At proper conditions, 
for example near the S-transition point, the potential wells can be very deep enabling the 
spontaneous generation of the vortex-antivortex pairs at the edges between potential val- 
leys and hills. The vortices and antivortices will screen these deep wells and hills similarly 
to the screening in the plasma. The difference is that, in contrast to plasma, the screen- 
ing "charges" do not exist without external potential. In such a flattened self-consistent 
potential relief the vortices have percolated infinite trajectories passing through the saddle 
points The drift motion of the delocalized vortices and antivortices in the external field 
generates dissipation and transfer the S-film into the resistive state (see figure llOKight). 
Replacing slow varying logarithmic potential by a constant at distances less than A and zero 
at larger distances, Feldman et al. have found thermodynamic and transport characteristics 
of this system. Below we briefly outline their main results. For the sake of simplicity we 
replace this slow varying potential V{y) by a potential having a constant value within the 
single cell: Vq = 2eo at the distance r < A and zero at r > A, where eo = $o/(16vr^A), $0 
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Figure 12: The static resistance p of the film vs dimensionless temperature t = T/Tc at 
typical values of parameters. 



is the magnetic flux quantum. Considering the film as a set of almost unbound cells of the 
linear size A we arrive at the following Hamiltonian for such a cell: 

H = -U^^airii + e'^n^i +2eo'^ninj, (61) 

i i i>j 

where rij is integer vorticity on either a dot and or a site of the dual lattice (between the 
dots) which we conventionally associate with location of unbound vortices. cTj = ±1, where 
subscript i relates to the dot, describes the random sign of the dot magnetic moments. cTj = 
on the sites of dual lattice. The first term of the Hamiltonian ()61|) describes the binding 
energy of the vortex at the magnetic dot and U ~ eo^d/^o; with $d being the magnetic 
fiux through a single dot. The second term in the Hamiltonian is the sum of single vortex 
energies, e = eo ln(Aefr/a), where a is the period of the dot array, ^ is the superconducting 
coherence length. The third term mimics the intervortex interaction. Redefining the constant 
e, one can replace the last term of equation by eo(I]ni)^. The sign of the vorticity on a 
dot follows two possible ('up'- and 'down-') orientations of its magnetization. The vortices 
located between the dots (n^ on the dual lattice) are correlated on the scales of order A and 
form the above-mentioned irregular checker-board potential relief. 

To find the ground state, we consider a cell with large number of the dots of each sign 
~ (A/a)^ 3> 1. The energy (jUT|) is minimal when the "neutrality" condition Q = X^'^j = 
is satisfied. Indeed, if Q 7^ the interaction energy grows as Q^, whereas the first term 
of the Hamiltonian behaves as \Q\ and can not compensate the last one unless Q ~ 1- 
The neutrality constraint means that the unbound vortices screen almost completely the 
"charge" of those bound by dots, that is K ^ (A^+ — N^) ~ \^N± ~ A/a where K is the 
difference between the numbers of the positive and negative dots and N± are the numbers of 
the positive and negative vortices, respectively. Neglecting the total charge |Q| as compared 
with A/a, we minimize the energy (j61|) accounting for the neutrality constraint. At Q = 



27 



the Hamiltonian ()61|) can be written as the sum of one-vortex energies: 



H = ^Hi] Hi = -U aiUi + enl. 



(62) 



The minima for any Hi is achieved by choosing rii = ni, an integer closest to the mag- 
nitude (jjK = cTj f//(2e). The global minimum consistent with the neutrality is realized by 
values of n-j that differ from the local minima values not more than over ±1. Indeed, in the 
configuration with rii = n°, the total charge is | ~ k\ X^cTiI = i^K. Hence, if k <C A/a, 

then the change of the vorticity at a small part of sites by ±1 restores neutrality. To be more 
specific let us consider i^' > 0. Let n be the integer closest to k, and consider the case k < n. 
Then the minimal energy corresponds to a configuration with the vorticity rii = —n at each 
negative dot and with the vorticity n or n — 1 at positive dots. The neutrality constraint 
implies that the number of positive dots with the vorticity n — 1 is M = Kn. In the opposite 
case K > n the occupancies of all the positive dots are n; whereas, the occupancies of the 
negative dots are either n or n + 1. Note that in our model the unbound vortices are absent 
in the ground state unless k is an integer. Indeed, the transfer of a vortex from a dot with 
the occupancy n to a dual site changes the energy by AE = 2e{K — n + 2). Hence, the 
energy transfer is zero if and only if k is an integer, otherwise the energy change upon the 
vortex transfer is positive. At integer k, the number of the unbound vortices can vary from 
to Kn without change of energy. The ground state is degenerate at any non-integer k 
since, while the total number of the dots with the different vorticities are fixed, the vortex 
exchange between two dots with the vorticities n and nil does not change the total en- 
ergy. Thus, our model predicts a step-like dependence of dot occupancies on k at the zero 
temperature and peaks in the concentration of unbound vortices as shown in figure ^2 The 
data for finite temperature were calculated in the Ref. The dependencies of the unbound 
vortex concentration on k for several values of x = e/T are shown in figure ITTl Oscillations 
are well pronounced for x ^ 1 and are suppressed at small x (large temperatures). At low 
temperatures, a; ^ 1, the half- widths of the peaks in the density of the unbound vortices are 
Ak, ^ 1/x and the heights of peaks are ~ jn, where 7 = K/N. 

Vortex transport- At moderate external currents j the vortex transport and dissipation 
are controlled by unbound vortices. The typical energy barrier associated with the vortex 
motion is eg. The unbound vortex density is m ~ a~^7 ~ (aA)~^ and oscillates with n as 
it was shown above. The average distance between the unbound vortices is / ~ y/aX. The 
transport current exerts the Magnus (Lorentz) force Fm = j^o/c acting on a vortex. Since 
the condition T ^ eo is satisfied in the vortex state everywhere except for the regions too 
close to Tc, the vortex motion occurs via thermally activated jumps with the rate: 



where /i = (^^cr„)/(47re^) is the Bardeen-Stephen vortex mobility The induced electric 
field is accordingly 




(63) 



Ec = lB/c = m^Qul/c, 



(64) 
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The Ohmic losses per unbound vortex are Wc = jE^Xa = J^quI/c giving rise to the dc 
resistivity as 

P-^- = -^ = ^ [-^0 (T) /T] (65) 

Note the non-monotonic dependence of pdc on temperature T figure ^1 The density of the 
unbound vortices is the oscillating function of the fiux through a dot. The resistivity of such 
a system is determined by thermally activated jumps of vortices through the corners of the 
irregular checkerboard formed by the positive or negative unbound vortices and oscillates 
with (^d- These oscillations can be observed by additional deposition (or removal) of the 
magnetic material to the dots. 



3.3 Ferromagnet - Superconductor Bilayer 
3.3.1 Topological Instability in the FSB 

Let us consider a F/S bilayer with both layers infinite and homogeneous. An infinite magnetic 
film with ideal parallel surfaces and homogeneous magnetization generates no magnetic field 
outside. Indeed, it can be considered as a magnetic capacitor, the magnetic analog of an 
electric capacitor, and therefore its magnetic field confined inside. Thus, there is no direct 
interaction between the homogeneously magnetized F-layer and a homogeneous S-layer in 
the absence of currents in it. However, Lyuksyutov and Pokrovsky argued [SO] that such a 
system is unstable with respect to spontaneous formation of vortices in the S-layer. Below 
we reproduce these arguments. 

Assume the magnetic anisotropy to be sufficiently strong to keep magnetization perpen- 
dicular to the film (in the z-direction). As we have demonstrated above, the homogeneous 
F-film creates no magnetic field outside itself. However, if a Pearl vortex somehow appears 
in the superconducting film, it generates magnetic field interacting with the magnetization 
m per unit area of the F-film. At a proper circulation direction in the vortex and the rigid 
magnetization m this field decreases the total energy over the amount m / Bz{r)(Px = m$, 
where $ is the total fiux. We remind that each Pearl vortex carries the fiux equal to the 
famous fiux quantum $o = vr^c/e. The energy necessary to create the Pearl vortex in the 
isolated S-film is 6^,°) = eoln(A/0 |Z21, where eo = ^g/lGyr^A, A = Xl/d is the effective 
penetration depth[71j, A^, is the London penetration depth, and ^ is the coherence length. 
Thus, the total energy of a single vortex in the FSB is: 

= el°) - m<l>o, (66) 

and the FSB becomes unstable with respect to spontaneous formation vortices as soon as e„ 
turns negative. Note that close enough to the S-transition temperature Tg, is definitely 
negative since the S-electron density and, therefore, e^^-* is zero at T,. If m is so small 
that e^, > at T = 0, the instability exists in a temperature interval < T < T^, where 
is defined by equation e^jiTy) = 0. Otherwise instability persists till T = 0. 

A newly appearing vortex phase cannot consist of the vortices of one sign. More accurate 
statement is that any finite, independent on the size of the film Lf density of vortices is 
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DOMAIN WALL 



Figure 13: Magnetic domain wall and coupled arrays of superconducting vortices with 
opposite vorticity. Arrows show the direction of the supercurrent. 



energetically unfavorable in the thermodynamic limit Lj ^ oo. Indeed, any system with 
the non-zero average vortex density n„ generates a constant magnetic field Bz = n^$o along 
the z direction. The energy of this field for a large but finite film of the linear size Lf 
grows as Lj exceeding the gain in energy due to creation of vortices proportional to L'j 
in thermodynamic limit. Thus, paradoxically the vortices appear, but can not proliferate 
to a finite density. This is a manifestation of the long-range character of magnetic forces. 
The way from this controversy is similar to that in ferromagnet: the film should split in 
domains with alternating magnetization and vortex circulation directions. Note that these 
are combined topological defects: vortices in the S-layer and domain walls in the F-layer. 
They attract each other. The vortex density is higher near the domain walls. The described 
texture represents a new class of topological defects which does not appear in isolated S and 
F layers. We show below that if the domain linear size L is much greater than the effective 
penetration length A, the most favorable arrangement is the stripe domain structure (see 
(figure Uni)). The quantitative theory of this structure was given by Erdin et al. ^T] . 
The total energy of the bilayer can be represented by a sum: 



where Ugv is the sum of energies of single vortices; Uyy is the vortex- vortex interaction energy; 
U^rn is the energy of interaction between the vortices and magnetic field generated by domain 
walls; Umm is the self-interaction energy of magnetic layer; Udw is the linear tension energy 
of domain walls. We assume the 2d periodic domain structure consisting of two equivalent 
sublattices. The magnetization m^(r) and density of vortices ?T,(r) alternate when passing 
from one sublattice to another. Magnetization is supposed to have a constant absolute value: 
mz{v) = ms(r), where s(r) is the periodic step function equal to at one sublattice and 
-1 at the other one. We consider a dilute vortex system in which the vortex spacing is much 
larger than A. Then the single- vortex energy is: 



U 




(67) 




'0 



(68) 
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The vortex-vortex interaction energy is: 



Uy^ = -Jn{r)V{r-r')n{r')d^xd^x', (69) 

where V{r — r') is the pair interaction energy between vortices located at points r and r'. 
Its asymptotics at large distances | r — r' |^ A is V^(r — r') = $0/ (^^r^ | r — r' |) jHSI- This 
long-range interaction is induced by magnetic field generated by the Pearl vortices and their 
slowly decaying currents'^. The energy of vortex interaction with the magnetic field generated 
by the magnetic film looks as follows [H^ : 

Uvm = I V</^(r - r')n(r') • a(™)(r)rfWx' (70) 

Here ip{r — r') = arctan ^^fr is a phase shift created at a point r by a vortex centered at 
a point r' and aS-'^\v) is the value of the vector-potential induced by the F-film upon the 
S-film. This part of energy similarly to what we did for one vortex can be reduced to the 
renormalization of the single vortex energy with the final result already shown in equation 
(|68p . The magnetic self-interaction reads: 

U^m = -f / B^:^\v)s{v)(fx (71) 

Finally, the domain walls linear energy is Udw = ^dwLdw where e^^ is the linear tension of 
the domain wall and L^w is the total length of the domain walls. 

Erdin et a/.jSl] have compared energies of stripe, square and triangular domain wall 
lattices, and found that stripe structure has the lowest energy. Details of calculation can 
be found in (see correction in i33^. The equilibrium domain width and the equilibrium 
energy for the stripe structure are: 



where m = m — e^/^Q and C=0. 57721 is the Euler constant. The vortex density for the stripe 
domain case is: 4 ~ X 

''^''^ = ~<!>lLsSm{7rx/L,) ^^^^ 

Note a strong singularity of the vortex density near the domain walls. Our approximation is 
invalid at distances of the order of A, and the singularities must be smeared out in the band 
of the width A around the domain wall. 



•^From this long-range interaction of the Pearl vortices it is ready to derive that the energy of a system of 
vortices with the same circulation, located with the permanent density Uy on a film having the lateral size 
L, is proportional to nlL^ 
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The domains become infinitely wide at T = T, and at T = T^. If e^^ < 4m^, tlie 
continuous approximation becomes invalid (see sec. 3.2.3) and instead a discrete lattice of 
vortices must be considered. It is possible that the long nucleation time can interfere with 
the observation of described textures. We expect, however that the vortices that appear first 
will reduce the barriers for domain walls and, subsequently, expedite domain nucleation. 

Despite of theoretical simplicity the ideal bilayer is not easy to realize experimentally. 
The most popular material with the perpendicular to film magnetization is a multilayerer 
made from Co and Pt ultrathin films (see Sec 13. 3. 4]) . This material has very large coercive 
field and rather chaotic morphology. Therefore, the domain walls in such a multilayer are 
chaotic and almost unmovable at low temperatures (see Sec 13. 3. 4|) . We hope, however, that 
these experimental difficulties will be overcome and spontaneous vortex structures will be 
discovered before long. 

3.3.2 Superconducting transition temperature of the FSB 

The superconducting phase transition in ferromagnet-superconductor bilayer was studied 
by Pokrovsky and Wei |33j. They have demonstrated that in the FSB the transition pro- 
ceeds discontinuously as a result of competition between the stripe domain structure in a 
FM layer at suppressed superconductivity and the combined vortex-domain structure in 
the FSB. Spontaneous vortex-domain structures in the FSB tend to increase the transition 
temperature, whereas the effect of the FM self-interaction decreases it. The final shift of 
transition temperature Tc depends on several parameters characterizing the SC and FM 
films and varies typically between -0.03Tc and 0.03Tc . 

As it was discussed earlier, the homogeneous state of the FSB with the magnetization 
perpendicular to the layer is unstable with respect to formation of a stripe domain structure, 
in which both, the direction of the magnetization in the FM film and the circulation of the 
vortices in the SC film alternate together. The energy of the stripe structure per unit area 
U and the stripe equilibrium width Lg is given in equations [721 IZ21 To find the transition 
temperature, we combine the energy given by equation [TSl with the Ginzburg-Landau free 
energy. The total free energy per unit area reads: 



Here a and /3 are the Ginzburg-Landau parameters. We omit the gradient term in the 
Ginzburg-Landau equation since the gradients of the phase are included in the energy (fTH)) . 
whereas the gradients of the superconducting electrons density can be neglected everywhere 
beyond the vortex cores. Minimizing the total free energy Pokrovsky and Wei [32] have 
found that near Tc the FSB free energy can be represented as 



— 16m^ 



exp(^-r^ + C - 1) + nsds[a{T - T^) + -n,] . 



(75) 



F = U + Fgl 
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where is given by the equation: 



= Tc H — exp(-— - + C - 1) (77) 



4^2 



The SC phase is stable if its free energy equation EHl is less than the free energy of a single 
FM film with the stripe domain structure, which has the following form [OHinSI: 

- (78) 



where Lf is the stripe width of the single FM film. Near the SC transition point the 
temperature dependence of the variation of this magnetic energy is negligible. Hence, when 
T increases, the SC film transforms into a normal state at some temperature T* below T^. 
This is the first-order phase transition. At transition point both energies Fg and Fm are 
equal to each other. The shift of the transition temperature is determined by a following 
equation: 



ATc = T^-Tc = — exp(— ^ + C - 1) - 



(79) 



Two terms in equation [7^1 play opposite roles. The first one is due to the appearance of 
spontaneous vortices which lowers the free energy of the system and tends to increase the 
transition temperature. The second term is the contribution of the purely magnetic energy, 
which tends to decrease the transition temperature. The values of parameters entering 
equation [7^1 can be estimated as follows. The dimensionless Ginzburg-Landau parameter 
is a = 7.04Tc/eir, where ep is the Fermi energy. A typical value of a is about 10~^ for 
low-temperature superconductors. The second Ginzburg-Landau parameter is (3 = aTcjUf,^ 
where ng is the electron density. For estimates Pokrovsky and Wei |^ take Tc ~ 3K, 
rig ~ 10^^ cm'^. The magnetization per unit area m is the product of the magnetization 
per unit volume M and the thickness of the FM film rf^- For typical values of M ~ 10^ 
Oe and dm ~ lOnm. m ~ 10"'* Gs/cm^. In an ultrathin magnetic film the observed 
values of Lf vary in the range 1 to 100/im EZI- If I// ~ l/xm, ds = dm = lOnm, and 

exp(-gd^/4m2 + C -1) ^ 10-^ ATjTc 0.03. For L/ = lOO/im, 4 = 50nm, and and 

exp(-l,,/4m2 + C - 1) ^ 10-^, ATJT, ~ 0.02. 



3.3.3 Transport properties of the FSB 

The spontaneous domain structure violates initial rotational symmetry of the FSB. There- 
fore, it makes transport properties of the FSB anisotropic. Kayali and Pokrovsky have 
calculated the periodic pinning force in the stripe vortex structure resulting from a highly 
inhomogeneous distribution of the vortices and antivortices in the FSB. The transport prop- 
erties of the FSB are associated with the driving force acting on the vortex lattice from an 
external electric current. In the FSB the pinning force is due to the interaction of the domain 
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Figure 14: Schematic vortex distribution in the FSB. The sign ± refers to the vorticity of 
the trapped flux. 



walls with the vortices and antivortices and the vortex-vortex interaction Uyy Periodic pin- 
ning forces in the direction parallel to the stripes do not appear in continuously distributed 
vortices. In the work [HI] the discreteness effects were incorporated. Therefore, one need to 
modify the theory [SI] to incorporate the discreteness effects. 

Kayali and Pokrovsky (33] have showed that, in the absence of a driving force, the vortices 
and antivortices lines themselves up in straight chains and that the force between two chains 
of vortices falls off exponentially as a function of the distance separating the chains. They also 
argued that pinning force in the direction parallel to the domains drops faster in the vicinity 
of the superconducting transition temperature and vortex disappearance temperature Ty. 

In the presence of a permanent current there are three kinds of forces acting onto a 
vortex. They are i) The Magnus force proportional to the vector product of the current 
density and the velocity of the vortex; ii) The viscous force directed oppositely to the vortex 
velocity; iii) Periodic pinning force acting on a vortex from other vortices and domain walls. 
The pinning force have perpendicular and parallel to domain walls components. In the 
continuous limit the parallel component obviously vanishes. It means that it is exponentially 
small if the distances between vortices are much less than the domain wall width. The sum 
of all three forces must be zero. This equation determines the dynamics of the vortices. It 
was solved under a simplifying assumptions that vortices inside one domain move with the 
same velocity. The critical current have been calculated for for parallel and perpendicular 
orientation. Theory predicts a strong anisotropy of the critical current. The ratio of the 
parallel to perpendicular critical current is expected to be in the range 10^ 10^ close to the 
superconducting transition temperature Tg and to the vortex disappearance temperature Ty. 
The anisotropy decreases rapidly when the temperature goes from the ends of this interval 
reaching its minimum somewhere inside it. The anisotropy is associated with the fact that 
the motion of vortices is very different in this two cases. At perpendicular to the domains 
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Figure 15: Magnetic properties of the Co/Pt multilayer: (a) Hysteresis loop measured 
by magneto-optical Kerr effect with H perpendicular to the sample surface. MFM images 
(5x5 /im^) show that the domain structure of the sample consists of band domains after 
out-of-plane demagnetization (b), bubble domains in the s = 0.3 (c) and s = 0.93 (d) states. 
(From Lange et al. ^cond-mat/0310132j . 

direction of the permanent current all the vortices are involved by the friction force into 
a drift in the direction of the current, whereas the Magnus force induces the motion of 
vortices (antivortices) in neighboring domains in opposite directions, both perpendicular to 
the current. The motion of all vortices perpendicular to the domains captures domain walls, 
which also move in the same direction. This is a Goldstone mode, no perpendicular pinning 
force appears in this case. The periodic pinning in the parallel direction and together with 
it the perpendicular critical current is exponentially small. In the case of parallel current 
the viscous force involves all vortices into the parallel motion along the domain walls and in 
alternating motion perpendicularly to them. The domain walls remain unmoving and provide 
very strong periodic pinning force in the perpendicular direction. This anisotropic transport 
behavior could serve as a diagnostic tool to discover spontaneous topological structures in 
magnetic-superconducting systems. 

3.3.4 Experimental studies of the FSB 

In the preceding theoretical discussion we assumed that the magnetic film changes its magne- 
tization direction in a weak external field and achieves the equilibrium state. All experimental 
works have been done with the Co/Pt, Co/Pd multilayers, which have large coercive field 
and are virtually "frozen" at the experiment temperature. Lange et al. pHl 1^ FlOOj have 
studied phase diagram and pinning properties of such magnetically "frozen" FSB. In these 
works the average magnetization is characterized by the parameter s, the fraction of spins 
directed up. Magnetic domains in Co/Pd(Pt) multilayers look like meandering irregular 
bands at s = 0.5 (zero magnetization) (see figure HHb) and as "bubble" domains (see figure 
ITRi) with typical size 0.25/im-0.35/im near fully magnetized states (s = or s = 1). The 
stray field from domains is maximal at s = 0.5 and decreases the superconducting transition 
temperature Tc of the Pb film by 0.2K (see figure IT^ . The effective penetration depth is 
about 0.76/im at 6.9K. 

Close to s = or s = 1 Lange et al. [HHl IHni HOOj have observed behavior in the 
applied magnetic field which is similar to the array of magnetic dots with normal to the film 
magnetization (see Sec. I3.2.1|) . They have found asymmetry in the applied magnetic field 
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Figure 16: Dependence of the critical temperature at zero field Tc(H=0) on the parameter 
s. The minimum value of is observed for s = 0.5. (From Lange et al. |cond-mat/0310132| ) 



for Tc{H) dependence and for pinning properties. The bubble domains have a perpendicular 
magnetic moment. If the thickness and magnetization is sufficient, they can pin vortices 
which appear in the applied external magnetic field. In this respect they are similar to 
randomly distributed dots with normal magnetization. Thus, in the range of filling factor 
s ~ Oorl the critical current must be large enough. Contrary to this situation at s ~ 0.5 
the randomly bent band domains destroy a possible order of the vortex lattice and provide 
percolation "routes" for the vortex motion. Thus the pinning is weaker and corresponds 
either to smaller critical current or to a resistive state. This qualitative difference between 
magnetized and demagnetized state has been observed in the experiments by Lange et al(|981 
inni llOOj ). The above qualitative picture of vortex pinning is close to that developed by 
Lyuksyutov and Pokrovsky j^E] and by Feldman et al [211 fo'^ transport properties of 
the regular array of magnetic dots with the random normal magnetization (see Sec 13. 2. 3]) . 
In this model the demagnetized state of the dot array is associated with the vortex creep 
through the percolating network. The strongly magnetized state, on the contrary, provides 
more regular vortex structure and enhances pinning. 

3.3.5 Thick Films 

Above we have discussed the case when both magnetic and superconducting films are thin, 
namely, ds Xl and dm ^ Lf. In this subsection we briefiy discuss, following works by 
Sonin [IH], situation when both films are thick ds ^ Xl and dm ^ Lf. Below we neglect 
the domain wall width. Consider first the ferromagnetic film without superconductor. This 
problem has been solved exactly by Sonin |101j . Figure IT7k shows schematically magnetic 
field distribution around thick ferromagnetic film. The problem can be solved by calculating 
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Figure 17: Magnetic charges (+ and -) and magnetic flux (thin hues with arrows) in a ferro- 
magnetic film (FM) without (a) and with (b) a superconducting substrate (SC). The mag- 
netization vectors in domains are shown by thick arrows. (From Sonin _cond-mat /0102102J 



field from "magnetic charges" on the magnetic film surface 

The magnetic field, without a superconducting substrate, at the ferromagnetic film 
boundary ?/ = 0+ is given by 



HJx) = -AM In 



tan 



nx 



2L, 



^0) 



Hy{x) = =F27rA'fsign tan 



nx 



2L 



at ?/ ±0 . 



^1) 



The field pattern is periodic with the period 2Lf along the axis x. The magnetic charge on 
the film boundary y = is 



pM = —M5{y)sign tan 



2Z7 



(82) 



Sonin has argued jlH| that in the case of bulk superconductor and with additional re- 
quirement Al/Lj — > 0, the magnetic fiux from the magnetic film is practically expelled 
from superconductor and problem can be solved by using images of magnetic charges on the 
magnetic film surface as shown in figure [T7b . Sonin has calculated energy change due to 
presence of the superconducting substrate and concluded that the substrate increases the 
total magnetic energy by 1.5 times. The energy of the domain walls per unit length along 
the axis x is inversely proportional to domain width Lfg and the energy of the stray fields 
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is proportional to Lfg. The domain width Lf^ is determined by minimization of the total 
energy per unit length. The growth of the magnetic energy decreases the domain width L fg 
by \/L5 times. Relative correction to the energy for finite X^/Lfg is of the order of \lI Lfg 

4 Proximity Effects in Layered Ferromagnet - Super- 
conductor Systems 

4.1 Oscillations of the order parameter 

All oscillatory phenomena theoretically predicted and partly observed in the S/F layered 
systems are based on the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) effect first proposed for 
homogeneous systems with coexisting superconductivity and ferromagnetism [23 IM] • They 
predicted that the energy favorable superconducting order parameter in the presence of 
exchange field should oscillate in space. The physical picture of this oscillation is as follows. 
In a singlet Cooper pair the electron with the spin projection parallel to the exchange field 
acquires the energy — /i, whereas the electron with the antiparallel spin acquires the energy 
+h. Their Fermi momenta therefore split onto the value q = 2h/vp. The Cooper pair 
acquires such a momentum and therefore its wave function is modulated. The direction of 
this modulation vector in the bulk superconductor is arbitrary, but in the S/F bilayer the 
preferential direction of the modulation is determined by the normal to the interface (z-axis). 
There exist two kinds of Cooper pairs differing with the direction of the momentum of the 
electron whose spin is parallel to the exchange field. The interference of the wave functions 
for these two kinds of pairs leads to the standing wave: 

F{z)=Fo cos qz (83) 

A modification of this consideration for the case when the Cooper pair penetrates to a 
ferromagnet from a superconductor was proposed by Demler et a/. |lU2j . They argued that 
the energy of the singlet pair is bigger than the energy of 2 electrons in the bulk ferromagnet 
by the value 2h (the difference of exchange energy between spin up and spin down electrons). 
It can be compensated if the electrons slightly change their momentum so that the pair will 
acquire the same total momentum q = 2h/vF- The value Im = vp/h called the magnetic 
length is a natural length scale for the LOFF oscillations in a clean ferromagnet. Anyway, 
equation (j83|) shows that the sign of the order parameter changes in the ferromagnet. This 
oscillation leads to a series of interesting phenomena that will be listed here and considered 
in some details in next subsections. 

1. Periodic transitions from 0- to vr phase in the S/F/S Josephson junction when varying 
thickness dj of the ferromagnetic layer and temperature T. 

2. Oscillations of the critical current vs. df and T. 

3. Oscillations of the critical temperature vs. thickness of magnetic layer. 
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Figure 18: Four types of trajectories contributing (in the sense of Feynman's path integral) 
to the anomalous wave function of correlated quasiparticles in the ferromagnetic region. 
The solid lines correspond to electrons, the dashed lines — to holes; the arrows indicate the 
direction of the velocity. (From Fominov et al. cond-mat / 0202280| ) 



The penetration of the magnetized electrons into superconductors strongly suppresses the 
superconductivity. This obvious effect is accompanied with the appearance of magnetization 
in the superconductor. It penetrates on the depth of the coherence length and is directed 
opposite to magnetization of the F-layer. Another important effect which does not have 
oscillatory character and will be considered later is the preferential antiparallel orientation 
of the two F-layers in the S/F/S trilayer. 

The described simple physical picture can be also treated in terms of the Andreev reflec- 
tion at the boundaries [il03j . long known to form the in-gap bound states [HS]; |104j . Due to 
the exchange field the phases of Andreev reflection in the S/F/S junction are different than in 
junctions S/I/S or S/N/S (with non-magnetic normal metal N). Indeed, let consider a point 
P inside the F-layer at a distance z from one of the interfaces jl()5j . The pair of electrons 
emitted from this point at the angle ib^, ±(7r — 6') to the 2-axis will be reflected as a hole along 
the same lines and returns to the same point f figure ITHj) . The interference of the Feynman 
amplitudes for these 4 trajectories creates an oscillating wave function of the Cooper pair. 
The main contribution to the total wave function arises from a small vicinity of = 0. Tak- 
ing only this direction, we find for the phases: 5*1 = —5*2 = —qz; = —S4 = —q{2df — z). 
Summing up all Feynman's amplitudes e*'^'=; k = 1...4, we find the spatial dependence of the 
order parameter: 

F (X cos qdf cos q{df — z) (84) 

At the interface F oc {cos qdj^. It oscillates as a function of magnetic layer thickness with 
the period Adf = Tc/q = 27TVF/h and decays due to the interference of trajectories with 
different 6. 

In a real experimental setup the LOFF oscillations are strongly suppressed by the elastic 
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impurity scattering. The trajectories are diffusive random paths and simple geometrical 
picture is not more valid. However, as long as the exchange field h exceeds or is of the 
same order of magnitude as the scattering rate in the ferro magnet 1/tj, the oscillations do 
not disappear completely. Unfortunately, the experiments with strong magnets possessing 
large exchange fields are not reliable since the period of oscillations goes to the atomic 
scale. Two layers with different thickness when they are so thin can have different structural 
and electronic properties. In this situation it is very difficult to ascribe unambiguously the 
oscillations of properties to quantum interference. 

4.2 Non-monotonic behavior of the transition temperature. 

This effect was first predicted by Radovic et al. t63.^. Its reason is the LOFF oscillations 
described in subsection 14.11 If the transparency of the S/F interface is low, one can expect 
that the order parameter in the superconductor is not strongly influenced by the ferro- 
magnet. On the other hand, the condensate wave function at the interface in the F-layer 
F oc (cos2(i//Am)^ becomes zero at df = nXm{n + l/2)/2 {n is an integer). At this values 
of thickness the discontinuity of the order parameter at the boundary and together with it 
the current of Cooper pairs into the ferromagnet has a maximum. Therefore one can expect 
that the transition temperature is minimal jlU6j . Experimental attempts to observe this 
effect were made many times on the S/F multilayers Nb/Gd |107j . Nb/Fe |lU8j . V/V-Fe 
[109 , V/Fe |110j . More references and details about these experiments and their theoretical 
description can be found in the cited reviews jHOllHI- Unfortunately, in these experiments 
the magnetic component was a strong ferromagnet and, therefore, they faced all the difficul- 
ties mentioned in subsection 13.3.11 the F-layer must be too thin and its variation produce 
uncontrollable changes in the sample, the influence of the growth defects is too strong. Be- 
sides, in multilayers the reason of the non- monotonous dependence of Tc on df may be the 
— vr transition. Therefore, the reliable experiment should be performed with a bilayer 
possessing a sufficiently thick F-layer. Such experiments were performed recently |11HI112] . 
The idea was to use a weak ferromagnet (the dilute ferromagnetic alloy Cu-Ni) with rather 
small exchange field h to increase the magnetic length Im = \J Df/h. They performed the 
experiments with S/F bilayers to be sure that the non-monotonic behavior is not originated 
from the — vr-transition. In these experiments the transparency of the interface was not 
too small or too large, the exchange field was of the same order as the temperature and the 
thickness of the F-layer was of the same order of magnitude as magnetic length. Therefore, 
for the quantitative description of the experiment theory should not be restricted by limiting 
cases only. Such a theory was developed by Fominov et a/. |105j . In the pioneering work by 
Radovic et al. jHHl the exchange field was assumed to be very strong. 

As always when it goes about critical temperature, the energy gap and anomalous Green 
function F are infinitely small. Therefore one needs to solve linearized equations of super- 
conductivity. The approach by Fominov et al. is based on solution of the linearized Usadel 
equation and is valid in the diffusion limit r^Tg ^ l;TfTc ^ ^ 1. Namely this sit- 
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Figure 19: FS bilayer. The F and S layers occupy the regions —dj < z < and < z < d, 
respectively. 



nation was realized in the cited experiments |11H 1111^ . The work by Fominov et al. jlU5j 
covers numerous works by their predecessors (HHl I113| I114t I1U2[ llU6j clarifying and improving 
their methods. Therefore in the presentation of this subsection we follow presumably the 
cited work |lU5j and briefly describe specific results of other works. 

The starting point is the linearized Usadel equations for singlet pairing for anomalous 
Green functions Fg in the superconductor and Ff in the ferromagnet: 

d'^F 

Ds^-\uJn\Fs + ^ = ^]Q<z<ds. (85) 

Df-Q^ - i\uJn\ + thsgnuJn)Ff = 0;-df<z< 0. (86) 

Thus, we accept a simplified model in which A = in the ferromagnetic layer and /i = in 
the superconducting one. The geometry is schematically shown in figure ^1 
Equations ()86l85p must be complemented with the self-consistency equation: 



A(r)ln^ = vrT5:f^-F. 

n, \ \^n\ 



:uJn,r)], (87) 



where T^s is the bulk SC transition temperature, and with linearized boundary conditions 
at the interface: 

dFs dFf , , 

dFf 

Aaj^ = G,iF.M-Ffm, (89) 

where agj is the conductivity of the superconducting (ferromagnetic) layer in the normal 
state; Gb is the conductance of the interface and A is the area of the interface. We assume 
that the normal derivative of the anomalous Green function is equal to zero at the interface 
with the vacuum: 

dFf dF 

The condition of solvability of linear equations ()85p8(if87j) with the boundary conditions 
fl88l89l90|) determines the value of transition temperature Tc for the F/S bilayer. 
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The solution F 



in the F-layer satisfying the boundary condition ()90|) reads: 



Ff{u)n, z) = C{un) cosh[A;/„(2; + df)]; kjn 



\ 



ujnl + ihsgnUn 



(91) 



where C{uJn) is the integration constant to be determined from the matching condition at 
the F/S interface z = 0. From the two boundary conditions at z = ()88|89|1 it is possible to 
eliminate Ff and dFj/dz and reduce the problem to finding the function Fg from equation 
(jH^ and the effective boundary condition at z = 0: 



is 



dF. 



7 



dz 7fe + Bf{ujn^ 



-Fs-, 



(92) 



where = Dg/ i2'KTcs]l = af/os]^b = iAaf)/{^sGb) and Bf{ujn) = (A;/„^s tanh(A;/„(i/)) . 
Since kfn is complex the parameter Bf{u!n) and consequently the function Fs is complex. The 
coefficients of the Usadel equation (|H3j) are real. Therefore, it is possible to solve it for the 
real part of the function Fg traditionally denoted as F+(c<j„, z) = ^ {Fgi^Un, z) + Fs{—uJn, z)). 
The boundary condition for this function reads: 



is 



dFl 
dz 



7(7fe + ^B 



f) 



(93) 



where 



- ksnds ta,nh.{k sndg) and kgn = y^\- To derive this boundary condition we 
accept the function A{z) to be real (it will be justified later). Then the imaginary part of 
the anomalous Green function F~(u;„, z) obeys the homogeneous linear differential equation; 



dz^ 



k^ F- 

'^sn-^ s 



and the boundary condition 



dF7 

dz 



at z = do. Its solution is F, 



J, Z) = E{uJn) COsh[ksniz- 

dg)]. At the interface z = its derivative ^^|^=o is equal to —ksntanh{ksnds)F~{LUn, z = 0). 
Eliminating F~ and its derivative from real and imaginary parts of the boundary condition 
fl92|l . we arrive at the boundary condition (jHS))- Note that only F^ participates in the self- 
consistence equation (jHTj) . This fact serves as justification of our assumption on reality of 
the order parameter A(z). 

Simple analytic solutions of the problem are available for different limiting cases. Though 
these cases are unrealistic at the current state of experimental art, they help to understand 
the properties of the solutions and how do they change when parameters vary. Let us consider 
the case of very thin S-layer dg -C is- In this case the order parameter A is almost a constant. 
The solution of equation (|H^ in such a situation is F^{uJn, z) = — + «n cosh ksn{z — dg), 
where a„ is an integration constant. From the boundary condition ()93|) we find: 



2AW{Un) 



\uJn\ {Agn + W{uJn))'' 



(94) 
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Figure 20: Theoretical fit to the experimental data. 
.cond-mat/0202280j . 
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(From Fominov et al. 



where the coefficients Asn are the same as in equation ()93|1 . We assume that ksnds ^ 1- 
Then Asn ~ ^l^^sds = |^(^ + 1/2)- The function F^{un,z) almost does not depend on z. 
The self-consistence equation reads: 

ln5^ = 2E V (95) 



The summation can be performed explicitly in terms of digamma-functions (F) 



1 -^cs 




Possible oscillations are associated with the coefficients Bf. If the magnetic length A 



(96) 



Df/h is much less than df, then Bf ^ sj ^Jr exp(— — In the opposite limiting 
case Xm ^ df there are no oscillations of the transition temperature. Note that ln(Tcs/Tc) 
can be rather large ~ C,s/ds, i.e the transition temperature in the F/S bilayer with very thin 
S-layer can be exponentially suppressed. This tendency is reduced if the resistance of the 
interface is large (7^ ^ 1). 

In a more realistic situation considered in the work |lU5j neither of parameters dg/ ^s, df / Xmi 'Jy 7fe 
is very small or very large and an exact method of solution should be elaborated. The au- 
thors propose to separate explicitly the oscillating part of the functions F^^Un, z) and A(z) 
and the reminders: 

= cos ad + ^ cosh a d ^^^^ 

^ ^ cosq{z-ds) ^ coshg^(z-4) ,gg. 

cosqds " coshg^^ 
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where the wave-vectors q and as well as the coefficients of the expansion must be found 
from the boundary conditions and self-consistence equation. Equation results in rela- 
tions between coefficients of the expansion: 

fn "i i ; f;; '^i fnm "i i (99) 
\uJn\ + DsQ^ \uJn\ - D^Q^ 

Substituting the values of coefficients fn, fnm from equation to the boundary condition 
fl9H|l ■ we find an infinite system of homogeneous linear equations for coefficients 6 and 6m: 

q tan qds -W{un) , qmtanhqmds - W{uJn) 
^ L; I I n ^2 + fr^2 (100) 

\^n\ + ^sQ m=l \^n\ — J^sQm 

Equating the determinant of this system T> to zero, we find a relation between q and qm- It 
is worthwhile to mention a popular approximation adopted by several theorists [HHl 110^1 llU6j 
the so-called single-mode approximation. In our terms it means that all coefficients Sm,Tn = 
1,2... are zero and only the coefficient 6 survives. The system (jlUUp implies that it is only 
possible when the coefficients W^Un) do not depend on their argument Un- It happens indeed 
in the limit dg/^s 1 and h^T. For a more realistic regime the equation V = must be 
solved numerically together with the self-consistence condition, which turns into a system of 
equations: 

T 
In^ 



c 



c 




In^ = F(--^:l^) -F(-) (101) 



These systems were truncated and solved with all data extracted from the experimental 
setup used by Ryazanov et al. |112j . The only 2 ffiting parameters were h = ISOi^' and 
76 = 0.3. 

Figure 1201 demonstrates rather good agreement between theory and experiment. Various 
types of the curves Tc{df) are shown in figure 1^ Note that the minimum on these curves 
eventually turns into a plateau at = 0, the reentrant phase transition into the super- 
conducting state. Some of the curves have a well-pronounced discontinuity, which can be 
treated as the ffist order phase transition. The possibility of the ffist order transition to 
superconducting state in the F/S bilayer was ffist indicated by Radovic et al. 



4.3 Josephson effect in S/F/S junctions 

As we already mentioned the exchange field produces oscillations of the order parameter 
inside the F-layer. This effect in turn can change the sign of the Josephson current in 
the S/F/S junction compared to the standard S/I/S or S/N/S junctions. As a result the 
relative phase of the S-layers in the ground state is equal to vr (the so-called vr-junction). 
In the closed superconducting loop with such a junction spontaneous magnetic flux and 
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Figure 21: Theoretical fit to the experimental data. (From Fominov et al. 
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spontaneous current appear in the ground state. These phenomena were first predicted by 
Bulaevsky et al. for vr-junction independently on the way of its realization. Buzdin et 
al. jSni have first argued that such a situation can be realized in the S/F/S junction at a 
proper choice of its length. Ryazanov et al. j^Zl HH] have realized such a situation employing 
the weak ferromagnet Cua;Nii_a; as a ferromagnetic layer. A similar approach was used by 
Kontos et al. who used a diluted alloy PdNi. The weakness of exchange field allowed 
them to drive the oscillations and in particular the O-tt transition by the temperature at a 
fixed magnetic field. The success of this experiment have generated an extended literature. 
The theoretical and experimental study of this and related phenomena still are active. In 
what follows we present a brief description of relevant theoretical ideas and the experiments. 



4.3.1 Simplified approach and experiment 

Here we present a simplified picture of the S/F/S junction based on the following assump- 
tions: 

i) The transparency of the S/F interfaces is small. Therefore the anomalous Green function 
in the F-layer is small and it is possible to use the linearized Usadel equation. 

ii) The energy gap A inside each of the S- layers is constant and equal to Aqc''^*'^/^ (the sign 
— relates to the left S-layer, + to the right one). 

iii) A = in the F-layer and /i = in the S-layers. 

The geometry of the system is shown in figure ()19|). From the second assumption it follows 
that the anomalous Green function F is also constant within each of S-layers: F = —. — = — =. 

The linearized Usadel equation in the F-layer ()86|) has a following general solution: 



F(cu„, z) = «,e'^/"^ + /5„e-'=^"^ (102) 

where 
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(compare equation ()91|)). The boundary condition at the two interfaces follows from the 
second boundary condition of the previous section (j89|) in which Ff is neglected: 

dFf , , 

is^ = (104) 

The coefficients and /?„ are completely determined by the boundary conditions ()l()4j) : 

cos(^±%^) , , 

smh{Kfndf) 

where Qn = ^? , ■ Equation (IH^ for the electric current must be slightly modified 

to incorporate the exchange field h: 

j = ze7rTiV(0)D Y.(FdF - FdF), (107) 

n 

where F{un, z) = F*{—uJn, z). Note that at this transformation the wave vectors kfn remain 
invariant. After substitution of the solution ()102|1 we find that j = jc sin ip with the following 
expression for the critical current |114[ I115[ : 



Jc D 



((^n + Al)kfndfsinh{kfndf] 

uj„>0 



(108) 



where Rn is the normal resistance of the ferromagnetic layer and Tf, is the dimensionless pa- 
rameter characterizing the ratio of the interface resistance to that of the F-layer. Kupriyanov 
and Lukichev [HH] have found the relationship between and the the barrier transmission 
coefficient Dh{6) [6 is the angle between the electron velocity and the normal to the interface): 

As we explained earlier, the oscillations appear since kjn are complex values. U h ^ 2'kT 



LUn, then kfn ~ (1 + 2d7 oscillations are driven only by the thickness. It was very 
important to use a weak ferromagnet with exchange field h comparable to ttT. Then the 
temperature also drives the oscillations. In the Cu-Ni alloys used in the experiment jSZ] the 
Curie point was between 20 and 50K. Nevertheless, the ratio h/nT was in the range of 
10 even for the lowest T^. In this situation kfn does not depend on n for the large number 
of terms in the sum (|108|) . This is the reason why the sum in total is periodic function of 



df with the period Am = '^y'^Dj /h. The dependence on temperature is generally weak. 
However, if the thickness is close to the value at which jc turns into zero at T = 0, the 
variation of temperature can change the sign of jc- 
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Figure 22: (Upper) Schematic cross-section of the sample. (Lower) Left : critical current Ic 
as function of temperature for Cuo.48Nio.52 junctions with differ- ent F-layer thicknesses be- 
tween 23 nm and 27 nm as indicated. Right : model calculations of the tem- perature depen- 
dence of the critical current in an SFS junction. (From Ryazanov et al. _cond-mat / 0008364J 



In figure ()22b ) theoretical curves jc{T) from cited work |lU5j are compared with the ex- 
perimental data by Ryazanov et al. jSZl ESI- The curves are plots of the modulus of jc vs 
T. Therefore, the change of sign of jc is seen as a cusp on such a curve. At temperature of 
the cusp the transition from 0- to vr-state of the junction proceeds. The change of sign is 
clearly seen on the curve corresponding to df = 2'lnm. The experimental S/F/S junction 
is schematically shown in figure ()22b). The details of the experiment are described in orig- 
inal paper [HZ] and in reviews [7^1 110] • Not less impressive agreement between theory and 
experiment is reached by Kontos et al. jHS] (theory was given by T. Kontos) (see figure 1^ . 

Very good agreement with the same experiment was reached in a recent theoretical work 
by Buzdin and Baladie |116j who solved the Eilenberger equation. 

Zyuzin et al. |117j have found that in a dirty sample the amplitude of the Josephson current 
jc is a random value with an indefinite sign. They estimated the average square fluctuations 
of this amplitude for the interval of the F-layer thickness < df < \JdJt as: 

where A is the area of interface and g is its conductance per unit area. The fluctuations are 
significant when df becomes smaller than the diffusive thermal length JD/T. 
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Figure 23: Josephson coupling as a function of thickness of the PdNi layer (full circles). The 
critical current cancels out at lii? ~ 65 A indicating the transition from "0" to "vr" -coupling. 
The full line is the best fit obtained from the theory as described in the text. Insert shows 
typical I-V characteristics of two junctions with (full circles), and without (empty circles) 
PdNi layer. (From Kontos et al. cond-mat /020 1 1 04J . 
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4.3.2 Josephson effect in a clean system 

In a recent work by Radovic et al |118j considered the same effect in a clean S/F/S trilayer. A 
similar, but somewhat different in details approach was developed by Halterman and Olives 
|119j . The motivation for this consideration is the simplicity of the model and very clear 
representation of the solution. Though in the existing experimental systems the oscillations 
are not disguised by impurity scattering, it is useful to have an idea what maximal effect 
could be reached and what role plays the finite transparency of the interface. The authors 
employed the simplest version of theory, Bogolyubov-DeGennes equations: 

(111) 





where a means —a and the effective Hamiltonian reads: 

_ ( H,{r) - ah{v) A(r) \ 
^-1, A*(r) -ifo + aMr)j' ^^^^^ 

Ho{v) = -^V'~^^ + W{v) (113) 

In the last equation /x is the chemical potential and W{v) is the barrier potential: 

W{Y) = W[5{z + d/2)+5{z-d/2)]. (114) 

The assumption about exchange field h{Y) and the order parameter A(r) are the same as 
in the previous subsubsection. We additionally assume that the left and right S-layers are 
identical and semi- infinite extending from z = —oo to z = —d/2 and from z = d/2 to 
z = oo. Due to translational invariance in the {x, ?/)-plane the dependence of the solution 
on the lateral coordinates is a plane wave: 

= e*''ii"^(2) (115) 

There are 8 fundamental solutions of these equations corresponding to the injection of the 
quasiparticle or quasihole from the left or from the right with spin up or down. We will 
write explicitly one of them "^ilz), corresponding to the injection of the quasiparticle from the 
right. In the superconducting area z < —d/2 we will see the incident quasiparticle wave with 
the coefficient 1 and the normal wave vector fc"*", the reflected quasiparticle with the reflection 
coefficient and the normal wave vector — the reflected quasihole (Andreev reflection) 
with the reflection coefficient ai and the wave vector k_, where {k^f = ^{Ep ± 0^ - ^, 




Ep is the Fermi energy and ^ = J E"^ — |Ap. Thus the solution "^liz) at z < —d/2 reads: 
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where u and v are the bulk Bogolyubov-Valatin coefficients: u = ^J{1 + ^/E)/2; v = sJ{l-i/E)/2. 
In the F-laycr —d/2 < z < d/2 there appear transmitted and reflected electron and trans- 
mitted and reflected hole. Since according to our assumption A = in the F-layer, there is 
no mixing of the electron and hole. With this explanation we can write directly the solution 
'^i{z) in the F-layer: 



= (Cie^^^^ + Cse^^^^) ' + (C3e^^-"^ + W'i'^') [ ] , (117) 





where — ^J^{E{p + ah± E) — k|. Finally in the right S-layer z > d/2 only the trans- 
mitted quasiparticle and quasihole propagate: 

- c,e-^^ ( ) + d,e--- ( ) (ns) 

The value of all coefficients can be established by matching of solutions at the interfaces: 

*(±--0) = *(±-+0); _U,,„-_U,_„ = ^* (119) 

Other fundamental solutions can be found by symmetry relations: 

a2{(p) = ai{-ip);a3 = 04 = aiib^ = 61; 64 = 62, (120) 

where index 2 relates to the hole incident from the left, indices 3,4 relate to the electron and 
hole incident from the right. Each mode generates the current independently on others. The 
critical current reads: 

eAT ^ K + K ( ain a2n\ 

Here all the values with the index n mean functions of energy E denoted by the same symbols 
in which E is substituted by for example ^„ = iyJuj^^-\-~A?. We will not demonstrate 
here straightforward, but somewhat cumbersome calculations and transit to conclusions. 
The critical current displays oscillations originated from two different types of the bound 
states. One of them appears if the barrier transmission coefficient is small. This is the 
geometrical resonance. The superconductivity is irrelevant for it. Another one appears even 
in the case of ideal transmission: this is the resonance due to the Andreev reflection. When 
the transmission coefficient is not small and not close to 1, it is not easy to separate these two 
type of resonances and the oscillations picture becomes rather chaotic. The LOFF oscillations 
are better seen when transmission coefficient is close to 1 since geometrical resonances do 
not interfere. Varying the thickness, one observes periodic transitions from to 7r-state 
with the period equal to A//2 = 2TrvF/h. The lowest value of d at which — 7r-transition 
takes place is approximately A//4. The temperature changes this picture only slightly, but 
near the thickness corresponding the — 7r-transition the non-monotonic behavior of jc vs. 



50 



temperature including temperature driven — vr-transition can be found. 
An intermediate case between the diffusion and clean limits was considered by Bergeret et al. 
|12Uj . They assumed that the F-layer is so clean that hrf 3> 1, whereas TcTg -C 1. Therefore 
Usadel equation is not valid for the F-layer and they solved the Eilenberger equation. They 
have found that the superconducting condensate oscillates as function of the thickness with 
period \f and penetrates into the F-layer over the depth equal to the electron mean free 
path If. The period of oscillations of the critical current is A//2. No qualitative differences 
with considered cases appear unless the magnetization is inhomogeneous. Even very small 
inhomogeneity can completely suppress the — vr-transitions. This is a consequence of the 
generation of the triplet pairing, which will be considered later. 



4.3.3 Half-integer Shapiro steps at the — vr transition 

Recently Sellier et al. |121j have reported the observation of the Shapiro steps at the voltage 
equal to half- integer of the standard values Vn = nhuj/2e, where oj is the frequency of 
the applied ac current. Let us remind that the standard (integer) Shapiro steps appear 
as a consequence of the resonance between the external ac field and the time-dependent 
Josephson energy Ej = —^cos(p{t) where the phase is proportional to time due to external 
permanent voltage through the contact: ip{t) = 2eVt/h. Just in the — tt transition point 
jc turns into zero. Then the next term in the Fourier-expansion of the Josephson energy 
proportional to cos(2(y9) dominates. That means that the Josephson current is proportional 
to sm{2ip). Such a term leads to the Shapiro steps not only at integer, but also at half-integer 
values since the resonance now happens at {4eV/h) = uj. Normally the term with sin(2(y9) is 
so small that it was always assumed to vanish completely. The resonance hf method used 
by the authors had sufficient sensitivity to discover this term. 

The authors prepared the Nb/Cu52Ni48/Nb junction by the photolitography method. Curie 
temperature of the F-layer is 20K. The 2 samples they used had the thicknesses 17 and 
19 nm. The — vr transition was driven by temperature. The transition temperature in 
the first and second sample were 1.12 and 5.36K, respectively. The external ac current had 
the frequency uj = 800 kHz and amplitude about 18 fiA. The voltage current curves for 
df = VJnm and temperatures close to 1.12 K are shown in figure (j23)- The fact that the 
half-integer steps disappear at very small deviation from the transition temperature proves 
convincingly that it is associated with the — vr transition. 



4.3.4 Spontaneous current and flux in a closed loop 



Bulaevsky et al. [55j argued that a closed loop containing the vr-junction may carry a 
spontaneous current and fiux in the ground state. Below we reproduce their arguments. The 
energy of the closed superconducting loop depends on the total fiux $ through the loop: 

£(*) = -Aj,eo=^ + (122) 

where (/? = Jc is the critical current and L is the inductance of the loop. The first 
term in equation ()122j) is the Josephson energy, the second is the energy of magnetic field. 



51 



10 20 30 40 50 60 70 

I { ) 



Figure 24: Shapiro steps in the vohage-current curve of a 17 nm thick junction with an 
excitation at 800 kHz (amphtude about 18 fiA). Half-integer steps (n=l/2 and n=3/2) 
appear at the O-vr crossover temperature T*. Curves at 1.10 and 1.07 K are shifted by 10 
and 20 /lA for clarity. (From Sellier et al. cond-mat / 0406236| ) . 

sorry, this figure is 
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Figure 25: Real (upper) and schematic (low) picture of the network of five SFS junctions 
Nh — Cuo,4QNiQ^4 — Nb {dp = 19 nm), which was used in the phase-sensitive experiment. 
(From From Ryazanov et al. |Cond-mat /0 103240, ) . 

The location of the energy minimum depends on the the parameter k = ^^^^ ^ . If k is 
positive, there is only one minimum at (f = 0. IfA;<— 1, the only minimum is located 
again at (f = 0. If— 1<A;<0, the minimum is located at the nonzero root of equation 
sm(f/(f = \k\; the value ip = corresponds to a maximum of energy. Thus, the spontaneous 
flux appears at sufficiently large inductance of the loop. It is possible to avoid this limitation 
measuring the dependence of the current inside the loop on the external flux through it |58j. 

They used triangular bridge array with vr-junctions in each shoulder (see Figs. 1^^ . Due 
to the central 7r-junction the phases of the current in two sub-loops of the bridge differ by 
TT. Therefore the critical current between the two contacts of the bridge is equal to zero in 
the absence of magnetic field. If the fiux inside the loop reaches half of fiux quantum, it 
compensates the indicated phase difference and the currents from both sub-loops are in the 
same phase. Thus, the shift of the current maximum from $ = to $ = $o/2 is the direct 
evidence of the — tt transition. Such experimental evidence was first obtained in the same 
work ^7\ . 

The graphs of the current vs magnetic field for two different temperatures (figure OBj) 
clearly demonstrates the shift of the current maximum from zero to non-zero magnetic field. 
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Figure 26: Magnetic field dependences of tlie critical transport current for the structure 
depicted in figure 1^ at temperature above (a) and below (b) Tcr- (From From Ryazanov et 
al. cond-mat/0103240). 
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Figure 27: (a) Temperature dependence of the critical transport current for the structure 
depicted in figure |2S1 in the absence of magnetic field; (b) temperature dependence (jump) 
of the position of the maximal peak on the curves Im{H), corresponding to the two limiting 
temperatures depicted in figure ESI (From From Ryazanov et al. cond-mat/0103240). 



The next graph figure 1771 shows the shift of the flux through the loop to 1/2 of the flux 
quantum at the temperature driven — vr transition. 

4.4 F/S/F junctions 

The trilayers F/N/F (N is normal non- magnetic metal) have attracted much attention start- 
ing from the discovery by Griinberg |122j of the Giant Magnetoresistance (GMR). The di- 
rection of magnetization of ferromagnetic layers in these systems may be either parallel or 
antiparallel in the ground state oscillating with the thickness of the normal layer on the scale 
of few nanometers. The mutual orientation can be changed from antiparallel to parallel by 
a rather weak magnetic field. Simultaneously the resistance changes by the relative value 
reaching 50%. This phenomenon has already obtained a technological application in the mag- 
netic transistors and valves used in computers |123j . A natural question is what happens if 
the central layer is superconducting: will it produce the spin-valve effect (a preferential mu- 
tual orientation of F-layers magnetization) and how does it depend on thicknesses of S and 
F-layers? This question was considered theoretically by several authors lUS} lU^l IH^l I124j . 
Recently the spin- valve effect was experimentally observed by Tagirov et al. |125j . 
Even without calculations it is clear that, independently on the thicknesses of S and F layers, 
the antiparallel orientation of magnetizations in F-layers has always lower energy than the 
parallel one. It happens because the exchange field always suppresses superconductivity. 
When the fields from different layers are parallel, they enhance this effect and increase the 
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energy, and vice versa. The effect strongly depends on the interfaces transparency. If it is 
very small, the effect is weak. In the case of almost ideally transparent interfaces the ma- 
jority electrons with the preferential spin orientation can not penetrate from the F-layer to 
the S-layer deeper than to the coherence length C,s- Therefore, it is reasonable to work with 
the S-layer whose thickness does not exceed ^s- The choice of the material and thickness of 
F-layers is dictated by the requirement that they could be reoriented by sufficiently weak 
magnetic field. Thus, the coercive force must be small enough. We refer the reader to the 
original works for quantitative details. 

An alternative approach is to study the thermodynamics of the F/S/F trilayer at a fixed mu- 
tual orientation of magnetic moments. Such a study was performed by Baladie and Buzdin 
[124 for the case of very thin superconducting layer dg -C C,s- They considered Fg almost 
as a constant, but incorporated small linear and quadratic deviations and solved the lin- 
earized Usadel equation as it was shown in subsection 14.21 to find the critical temperature 
vs. thickness of the ferromagnetic layers. They have found that at large '-fb (low interface 
transparency) the transition temperature monotonically decreases with df increasing from 
its value in the absence of the F-layers to some saturation value and there is no substantial 
difference between parallel and antiparallel orientations. At smaller values of the suppres- 
sion of Tc increases and at parallel orientation the reentrant transition occurs at df ~ ^f, 
but still the transition temperature saturates at large df. At 76 smaller than a critical value 
the transition temperature becomes zero at a finite thickness df for both parallel and an- 
tiparallel orientation. The authors also have found some evidences that at low 7?, the SC 
transition becomes discontinuous for the parallel orientation. This conclusion was confirmed 
by a recent theoretical study by Tollis |12(ij . who has proved that the SC transition for the 
antiparallel orientation is always of the second order, whereas for the parallel orientation it 
becomes of the first order for small |126j . Baladie and Buzdin jl24j have considered also 
the energy gap at low temperature. For the case of thick ferromagnetic layers df ^ ^f they 
have found that the energy gap is the monotonically decreasing function of the dimensionless 
collision frequency (t/Aq)"^, where Aq is the value of the energy gap in the absence of the 
ferromagnetic layers. It turns into zero at (tjAq)"^ = 0.25 for the parallel and 0.175 for the 
antiparallel orientation. 

4.5 Triplet pairing 

If the direction of the magnetization in F-layer is inhomogeneous due to a domain wall 
or artificially, the singlet Cooper pairs penetrating into the F- from S-layer will be partly 
transformed into the triplet pairs. This effect was first predicted by Kadigrobov et al. [HZ] 
and by Bergeret et al. 6F. The triplet pairs cannot penetrate to the superconductors over 
the length larger than magnetic length = \l Df/h (or vp/h for the clean ferromagnet), but 
in the ferromagnet they are neither exchange interaction nor the elastic scattering suppresses 
them. Therefore, they can penetrate over much longer distance = \J Df/T. Even if the 
triplet pairing is weak, it provides the long-range coupling between two superconducting 
layers in a S/F/S junction. Moreover, if the thickness df exceeds Im significantly, only 
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triplet pairs survive at distances much larger than completely changing the symmetry 
properties of the superconducting condensate. 

The exchange field rotating in the y — 2;-plane is naturally described by the operator in the 
spin space h = hi^a^ cos a + a2sma), where /i is a scalar function of coordinates, (J2 and 
(T3 are the Pauli matrices and the angle a is a function of coordinates. It is clear that the 
non-diagonal part of h fiips one of spins of the pair transforming the singlet into the triplet. 
It does not appear if the magnetization is coUinear (a = 0). To make things more explicit, 
let consider the Usadel equation in the F- layer, i.e. equation of the Section (j2)). First 
we simplify them by linearization, which is valid if either the transparency of the interface 
barrier is small |H3I- Then the condensate Green tensor / in F-layer is small. The linearized 
Usadel equation reads: 

- + fo{(3-3,/cosa + f3[(3-2,/]sina} =0, (123) 

where {A, B} means the anticommutator of operators A and B. If a = const, equations 
have an exponential solution / = e^^/o. The secular equation for k is: 



2\2 



ik' - klf + ^ 



0, (124) 



where fc^ = 2\u\/Df. Note that the secular equation does not depend on a. It is a conse- 
quence of rotational invariance of the exchange interaction. At a = the two-fold eigenvalue 
k"^ = k"^ corresponds to /i^2 (triplet pairing with projection ±1 onto the magnetic field). Since 
Un is proportional to T, these modes are long-range. Two other modes have wave vectors 
k = kh and k = fc^, where k^ = 2{\lj\ + ihsignco) / D f . They penetrate not deeper than on the 
magnetic length. These short-range modes are linear combinations of the singlet and triplet 
with spin projection zero, i.e. orthogonal to the magnetic field. 

Bergeret et al. considered two different geometries. In the first one [HHI they considered S/F 
bilayer. The angle a was a linear function of coordinate starting from at the S/F interface, 
reaching a value at the distance w from the interface and remaining constant at larger dis- 
tances. They have solved the linearized Usadel equation p23p with the boundary condition 
^f~dz' ~ IbFs proper at small transparency of the interface by a clever unitary transforma- 
tion / ^ lJ{z)f[ij{z)]-'^ with lJ{z) = exp{iQaiz/2) and Q = ^ = ^. This transformation 
turns the rotating magnetic field into the constant one, directed along 2;-axis, but differential 
term generates perturbations proportional to Q and Q^. By this trick the initial equations 
with the coordinate dependent h{z) is transformed into an ordinary differential equation 
with constant (operator) coefficients. The generation of the triplet component is weak if 7f, 
is large and it acquires an additional small factor if the ratio ^j/w is small {w mimics the 
domain wall width), but, as we have demonstrated, this component has a large penetration 
depth. Experimentally it could produce a strong enhancement of the F-layer conductivity. 
Such an enhancement was observed in the experiment by Petrashov et a/. |127j in 1999, two 
years afore the theoretical works. They studied an F/S bilayer made from 40nm thick Ni 
and 55nm thick Al films. The interface was about lOOxlOOnm^. The samples were prepared 
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Figure 28: 6-layer structure. 



by e-beam lithography. They measured the resistivity and the barrier resistance directly. 
They have found also the diffusion coefficients Dg = lOOcm^/s and Df = lOcm^/s, which 
we cite here to give an idea about the order of magnitudes. They have found a large drop 
of the resistance of the sample, which could not be explained by the existing singlet pairing 
mechanism. We are not aware about the detailed comparison of the theory jHHl and the 
experiment jl27j . One more evidence of long range penetration of the superconducting or- 
der parameter through the ferromagnet was reported in jl28j . The authors measured the 
resistance of O.bfim Ni loop connected with superconducting Al wire. They extracted the 
decay length for proximity effect in ferromagnet from differential resistance and concluded 
that it is much larger than it could be expected for singlet pairing. 

In their second work on the triplet pairing |^ the authors have proposed an interesting 
6-layer structure presented in figure |2H1 The assume that the magnetization in each layer is 
constant, but its direction is different in different layers. It is supposed to lay in the y — z- 
plane and thus it can characterized by one angle. Let this angle is —a in the layer Fi, in 
the layer F2 and ±a in the layer F3. They speak about the positive chirality if the sign is -|- 
and negative chirality if the sign is — . They prove that, if the thickness of F-layers is larger 
than Im, the superconducting layers Sa and Sb are connected by 0-junction if the chirality is 
positive and by vr-junction if the chirality is negative. This phenomenon is completely due to 
the triplet pairing since it dominates on this distance. Kulic and Kulic |129j considered two 
bulk magnetic superconductors with rotating magnetization separated by an insulating layer. 
They also have found that the sign of the Josephson current can be negative depending on 
the relative chirality. In this system singlet and triplet pairs coexist in the bulk, whereas in 
the system proposed by Bergeret et al. the triplet dominates. We will give a brief description 
how did they derive their results. They solved the Usadel equation in each layer separately 
(it can be done without linearization, since the coefficients of the differential equations are 
constant) and match these solutions using the Kupriyanov-Lukichev boundary conditions. 
The current density in the F2 layer can be calculated using the modified Eilenberger- Usadel 
expression: 



The maximal effect is reached when magnetic moment of the central layer is perpendicular 
to two others. 




(125) 
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5 Conclusions 



This short review shows that though the studies of Ferromagnet-Superconductor Hybrids 
are coming of age, we are at the beginning of interesting voyage into this emerging field. The 
most active development undoubtedly take place in the field of proximity based phenomena 
in layered fcrromagnet-superconductor systems. The strong point of this thrust is fruitful 
collaboration between experiment and theory. This progress was achieved due to a new 
idea due to Ryazanov and coworkers to use the weak ferromagnets in th experiment. This 
idea allowed to increase the thickness of ferromagnetic layers to a macrosopic scale and 
simultaneously allowing to drive the non-monotonous behavior of the Josephson current by 
temperature. On this way experimenters have reliably found several interesting phenomena 
predicted many years ago, as — 7r-transition and oscillations of critical temperature vs. the 
thickness of the F-layer and also some new phenomena as the valve effect in F/S/F junction 
and the Shapiro steps at half- integer frequencies. 

The experimental studies of ordering/transport in FSH have greatly benefited with intro- 
duction of imaging technique (SHPM,MFM) in the field. We expect that several experimental 
groups will get access to this technique in the near future which will result in more exciting 
experiments. The theoretical and experimental studies of ordering/transport in FSH have 
surprisingly little overlap, especially in comparison with studies of proximity based phenom- 
ena. The materials used in the experiment are far from being regular, whereas the theorist 
so far preferred simple problems with regular, homogeneous or periodical systems. Even the 
simplest idea about topological instability in the S/F-bilayer was not checked experimentally. 
It would be very instructive to find experimentally the phase diagram of a single magnetic 
dot using the SQUID magnetometer or the MFM. Finally the transport properties of the 
S/F- bilayer and the S-films supplied with regular or randomly magnetized arrays of F-dots 
should be measured. On the other hand the experiment dictates new problems for theory: 
a description of random set of strongly pinned domain walls, their magnetic field and its 
effect on the S-films. We think that both experimental and theoretical communities can find 
systems of common interests. Another possibility for interesting development in the FSH 
field we expect with introduction of new types of FSH, e.g. arrays of magnetic nanowires in 
alumina templates, covered with superconducting film. Such arrays provides alternative to 
magnetic dots source of alternating magnetic field of high strength and short scale variation. 
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